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rKEFACE 


The object of this book is to present a logical development 
of electromagnetic theory founded upon the principle of rela- 
tivity. So far as the author is aware, the universal procedure 
has been to base the electrodynamic equations on the experi- 
mental conclusions of Coulomb, Ampere, and Faraday, even 
books on the principle of relativity going no farther than to 
show that these equations are covariant for the Loren tz -Einstein 
transformation. As the dependence of electromagnetism on the 
relativity principle is far more intimate than is suggested by 
this covariance, it has seemed more logical to derive the electro- 
dynamic equations directly from this principle. 

The analysis necessary for the development of the theory has 
been much simplified by the use of Gibbs’ vector notation. 
While it is difficult for those familiar with the many conven- 
iences of this notation to understand why it has not come into 
universal use among physicists, the belief that some readers 
might not be conversant with the symbols employed has led to 
the presentation in the Introduction of those elements of vector 
analysis which are made use of farther on in the text. 

Chapter I contains a brief account of the principle of rela- 
tivity. In the second chapter the retarded equations of the 
field of a point charge are derived from this principle, and in 
Chapter III the simultaneous field of a moving charge is dis- 
cussed in some detail. In the next chapter the dynamical equa- 
tion of the electron is obtained, and in Chapter V the general 
field equations are derived. Chapter VI takes up the radiation 
of energy from electrons, and Chapters VII and VIII contain 
some applications of the electromagnetic equations to material 
media, chosen as much for their illustration of the theory as 
for their fundamental importance. Throughout, great pains 
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have been taken to distinguish between definitions and assump- 
tions, and to carry on the physical reasoning as rigorously as 
possible. It is hoped that the book may be found useful by 
those lecturers and students of electrodynamics who are looking 
for a logical rather than a historical account of the science. 
The subject matter covers topics appropriate for a one-year 
graduate course in electrodynamics and electromagnetic theory 
of light. 

The author wishes to acknowledge his debt to those great 
thinkers, Maxwell, Poynting, Gibbs, Lorentz, Larmor, and Ein- 
stein, and to express his appreciation of the inspiration and un- 
failing interest of his former teacher, Professor H. A. Bumstead. 
His thanks are due his colleague, Professor H. S. Uhler, for many 
suggestions tending toward greater clearness of exposition. 


Yale University 
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AN INTRODUCTION TO 
ELECTRODYNAMICS 

INTRODUCTION 

ELEMENTS OF VECTOR ANALYSIS 

Addition and multiplication. A vector is defined as a quantity 
which has both magnitude and direction. It will be designated 
by a letter in blackface type, its scalar magnitude being repre- 
sented by the same letter in italics. Geometrically, a vector may 
be represented by an arrow having the direction of the vector 
and a length proportional to its magnitude. The beginning of this 
representative straight line is known as its origin, and the end, 
%s its terminus . To add two vectors P and Q place the origin 
Df Q at the terminus of P. Then the line drawn from the 
origin of P to the terminus of Q is defined as the sum of P 
and Q. To subtract Q from P reverse the direction of Q and 
add. The components of a vector are any vectors whose sum is 
equal to the original vector. Although, strictly speaking, the 
components of a vector are themselves vectors, the term com- 
ponent will often be used to denote the magnitude alone in 
cases where the direction has already been specified. 

A vector is often determined by its components along three 
mutually perpendicular axes X \ Y, Z. These axes will always 
be taken so as to constitute a right-handed set ; that is, so that 
a right-handed screw parallel to the Z axis will advance along 
this axis when rotated from the X to the Y axis through the 
right angle between them. . Let i, j, k be unit vectors parallel 

l 
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respectively to the X, Y, Z axes. Then if the projections of P 
along these axes are denoted by /*, P u , /', 

P=# +%}+&, (1) 

and, obviously, 

P + Q = « + i + 05 + ft) i + ( P * + ft) k - ( 2 ) 

If two or more vectors are parallel to the same straight line, 
they are said to be collinear. If three or more vectors are par- 
allel to the same plane, they are said to be coplanar. 

Two vectors P and Q may be multiplied together in three 
different ways. The most general type of multiplication yields 
the undetermined product given by 

PQ-Wi+Wj+^ftik 

+ + + (3) 

+ P,QJd + P.Q#+P'QJ&. 

This product is neither vector nor scalar ; it is known as a dyad . 

The vector or cross product of two vectors is a vector perpen- 
dicular to their plane in the direction of advance of a right- 
handed screw when rotated from the first to the second of these 
vectors through the smaller angle between them. Its magnitude 
is equal to the product of the magnitudes of the two vectors by 
the sine of the angle between them. Therefore 

PxQ = — QxP. (4) 

Geometrically, this vector product has the magnitude of the 
parallelogram of which P and Q are the sides, and a direction 
at right angles to its surface. It follows from simple geometrical 
considerations that the distributive law holds for this product, 
thatis, (? -f Q) X R= PxR + QxR. (5) 

Therefore, inserting crosses between the vectors in each term 
of (3), 

Px Q= ( ZQ-P'Qy) i + (m ~ £ft) J + 05ft - W k - ( 6 ) 

The scalar or dot product of two vectors is a scalar equal in 
magnitude to the product of the magnitudes of the two vec- 
tors by the cosine of the angle between them. Obviously the 
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distributive law holds for this product. Therefore, inserting 
dots between the vectors in each term of (3), 

V'Q=P x Q x +PyQ»+ P zQz- ( 7 ) 

The triple scalar product 

(P x Q)-R 

evidently measures the volume of the parallelepiped of which 
P, Q, and R are the edges. Hence the position of cross and 
dot in this product is immaterial, and its sign is changed by 
interchanging the positions of two adjacent vectors. 

The triple vector product 

(P x Q) x R 

is obviously a vector in the plane of P and Q. From simple 
geometrical considerations it follows that 

(P x Q) x P = P 2 Q — P-QP| (8) 

and (P x Q) x Q = P-QQ - Q 2 P. (9) 

Now R may be written 

R = aP + JQ + c?(P x Q). 

Therefore 

(P x Q) x R = a(P x Q) x P + b(P * Q) x Q 

= OP 2 + iP-Q) Q - («P-Q + ^0 2 ) P 
= P-RQ - Q-RP. (10) 


This important expansion may be put in words as follows: 
Dot the exterior vector into the remoter vector inside the pa- 
rentheses to form the scalar coefficient for the nearer one, then 
dot the exterior vector into the nearer vector to form the scalar 
coefficient for the remoter one, and subtract this result from 
the first. 

The vector operator V (read del) is of great importance in 
mathematical physics. This quantity is defined as 




4 AN INTRODUCTION TO ELECTRODYNAMICS 


Let <f> be a scalar function of position in space. Then 

(11) 

is known as the gradient of <j>. It may easily be shown to rep- 
resent both in magnitude and direction the greatest (space) rate 
of increase of 4> at the point in question. 


Let V = Ki + VJ + V z k 

be a vector function of position in space. Then 


V-V = 


dx cy 



( 12 ) 


is known as the divergence of V. If V is the flux of a fluid 
per unit time per unit cross section, the divergence of V is the 
excess of flux out of a unit volume over that into this volume. 
If the fluid is incompressible, the divergence is obviously zero 
except at those points where sources or sinks are present. 

The vector 


VxV= 





(13) 


is known as the curl of V. If V specifies the linear velocities 
of the points of a rigid body, the curl is equal in magnitude and 
direction to twice the angular velocity of rotation. 

The following identities may easily be verified by expansion : 


V x V<£ == 0, (14) 

V-V x V== 0, (15) 

V x V x V = VV»V — V*VV. (16) 


Gauss 1 Theorem. In treating vector integrals volume, surface, 
and line elements will be denoted respectively by cZt, cfcr, and d \ . 
The direction of an element of a closed surface will be taken 
as that of the outward-drawn normal, and the direction of an 
element of a closed curve will be taken as that in which a right- 
handed screw passing through the surface bounded by the curve 
must rotate in order to advance toward the positive side of 
this surface. 
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Let V be a vector function of position in space, 
theorem states that 



Then Gauss’ 
( 17 ) 


where the surface integral is taken over the surface cr bounding 
the volume r. 

This theorem may be proved in the following way. In rectan- 
gular coordinates 



Let y , z and x 2 , y , z be the points of intersection of the 

surface bounding r with a line parallel to the X axis. Then 



Therefore 


dxdydz — J { V x (x 2 , y, z)~ V x (x v y, z)}dydz 

— Jv x dydz. 

jfv.v*=J (W*+ V y dzdx 4- V z dxdy) 

= jVcfrr. 


Stokes’ Theorem. If V is a vector function of position in 
space, Stokes’ theorem states that 



X V-rfcr 


=Jv.d\, 


( 18 ) 


where the line integral is taken over the curve \ bounding the 
surface cr. 

To, prove this theorem proceed as follows. In rectangular 
coordinates 
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Let x, y v z l and x, y 2 , z 2 be the points of intersection of the 
periphery of <r with a plane parallel to the YZ coordinate plane. 
Then, taking account of the signs of the differentials involved, 


J'(~dzdx — ^ dxdyj= J [{V x (x, y+dy, z + dz) 

- V x (x, y + dy , z')} dx + {V x ( x , y + dy, z) - V x (x, y, z)} dx] 


Therefore 



xV*c?«r 


= f {K(*, y? y? *,)}<& 

= f V x dx. 

= j(Vjx + V y dy + V,dz) 



Dyadics. A dyadic is a sum of a number of dyads. The first 
vector in each dyad is called the antecedent, and the second the 
consequent . Any dyadic may be reduced to the sum of three 
dyads. For if the dyadic v|/ is given by 

v|/ = al + bm + cn + do, (19) 

the vector o may be written 

o =/l + gtn + An, 

whence \|/ = (a -f/d) 1 4- (b -f #d) m -f (c 4- Ad) n. (20) 

Similarly, if either the antecedents or consequents of a dyadic 
are coplanar, the dyadic may be reduced to the sum of two 
dyads. Such a dyadic is said to be planar. If either anteced- 
ents or consequents are coll inear, the dyadic becomes a single 
dyad and is said to be linear . 

Consider the dyadic 

\|/ = al + bm + cn. 

If P is a vector, 

i|r»P = al*P +bm*P + cn»P 

is also a vector. Dotting a dyadic into a vector, then, gives rise 
to a vector having a new direction and magnitude. This new 



ELEMENTS OF VECTOR ANALYSIS 


7 


vector is a linear vector function of the original one. If a dyadic 
is planar, it will reduce to zero vectors having a certain direc- 
tion, and if it is linear, it will cause all vectors parallel to a 
certain plane to vanish. 

Obviously any dyadic may be written in the expanded form 
y\i = a u ii + a u ij + a ls ik 

+ a J i + tt Jj + a 2 J k ( 21 ) 

+ a n M + a M k l + 


It will now be shown that any dyadic may be put in such a 
form that its antecedents and consequents each constitute a 
right-handed set of mutually perpendicular vectors. Let a be 
a unit vector of variable direction extending from the origin. 

Then p = \|/»a 

describes a closed surface about the origin as a varies in direction. 
This surface may easily be shown to be an ellipsoid. Let i be 
the value of a for which p assumes its maximum value a. Now 
consider all values of a lying in the plane perpendicular to i. 
Let j be the value of a in this plane for which p assumes its 
greatest value b. Finally, let k be a unit vector perpendicular 
to i and j in the sense that will make i, j, k a right-handed set. 
Let c be the value of p when a equals k. Then, as the dyadic 
changes i, j, k into a, b, c, it may be written in the form 

\|r = ai + b j + ck. 

Now p = (ai + bj + ck)-a, 

d p = (ai + bj + ck)*da, 
p-dp = p-ai«da # + p*bj*da + p*ck*da. 


When a is parallel to i, p has its maximum value a, and 
therefore a-bj-tfa + a-ck-da = 0. 


If, moreover, da is perpendicular to j, a*c vanishes, and if da 
is perpendicular to k, a*b vanishes. Hence a is perpendicular 
to both b and c. 
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Now let a be restricted to the jk plane. Then 

= p*bj«c?a + P*ck«da. 

When a is parallel to j, p has its greatest value b, and therefore 
b»c vanishes. Therefore a, b, and c are mutually perpendicular. 
If they do not form a right-handed set, the direction of one of 
them may be reversed provided its sign is changed. Hence, if 
i x , ]\, kj constitute a right-handed set of mutually perpendicular 
unit vectors parallel respectively to a, b, c, the dyadic may be 
written v|/ = aiii + jjj + ^k. (22) 

If }\f = al + bm + cn, 

the conjugate of \| t is defined as 

i|/ c =la + mb + nc, 

and a dyadic is said to be self-conjugate if it is equal to its 
conjugate. Obviously, if \f/ is self-conjugate, 

\(/.P = P-v| r. (23) 

The idemf actor I is defined as 

I = ii + j j + kk 

= Mi + jii + kAv 

Evidently this dyadic is self-conjugate, and moreover 

I-P = P. (24) 

It will now be shown that any complete self-conjugate dyadic 
may be written in the form 

\|/ = aii + b)] -f c‘kk, (25) 

where i, j, k constitute a right-handed set of mutually perpen 
dicular unit vectors. 

It has already been shown that the dyadic may be written 
\|/= ai ± i + JjJ + cls^k. 

Now, as \|/ is self-conjugate, 

= aiii 4- ijj t + ckk v 
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Form the products 

+ * a jJi + <?KK> 

and >\f c *y\r = a 2 ii + b 2 jj + c 2 kk. 

These are equal, as the dyadic is self-conjugate. Therefore 
put a|/ 2 = \|/.\|/ c = 

Now consider the dyadic 

y\f 2 — a 2 l = (b 2 — a 2 ) jj*! -f (V — a 2 ) k x k r 
Obviously (\|/ 2 — a 2 T)'i x = 0. 

But, as \|/ a — a 2 I = (ft 2 — a 2 ) j j H- (c 2 — a 2 ) kk, 
it follows that (\|/ 2 — a 2 I)*i = 0. 

Hence, as (+■-<«) is planar, but not linear, i x must be 
parallel to i. Similarly, j x must be parallel to j, and k x to k. 



CHAPTER I 


THE PRINCIPLE OF RELATIVITY 

1. Motion. The concept of motion comprises two essential 
factors : a moving element , and a reference body relative to which 
the motion takes place. A grain of sand lying on the floor of a 
railway car is not in motion at all if the car itself is chosen as 
reference body, although it may be moving rapidly relative to 
the earth. If, however, the ear, the earth, and all other objects 
save the grain of sand are removed, the lack of a reference body 
makes it impossible to form a conception of motion. 

A moving element is characterized by a point — whether in 
a material body or not — which can be continuously identified. 
In the following discussion a point always will be understood 
to have this property. A reference body is essentially a group 
of points along the path of a moving element, together with a 
device for assigning numerical values to tlie intervals of time 
between coincidences of the moving element with successive 
points of the body, and to the distances between these points. 
For such characteristics are obviously necessary in order to 
make possible quantitative evaluation of the motion of the 
moving element. 

2. Reference system. A reference system is an assemblage of 
points filling all space. A device is provided for indicating 
time at these points in such a way as to assure synchronism 
according to some arbitrary standard, and for measuring distances 
between them. This device is subject to the following con- 
ditions, but otherwise it is quite arbitrary: 

(1) Two points which are in synchronism with a third are 
also in synchronism with each other. 

(2) The distance between two points is independent of the 
time at which it is measured. 


10 
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Thus a reference system serves as a reference body for any 
moving element. It must not, however, be imagined to offer 
any obstruction to the motion through it of such an element, 
or of another reference system. 

A material body of finite extent may be considered to consti- 
tute a reference system if the points of the body itself are sup- 
posed to have points outside associated with them in such a 
way that the whole assemblage possesses the properties described 
above. In order that the material part of such a system shall in 
no degree obstruct the motion through it of a moving element, 
those portions of it which would be in the way may be regarded 
as temporarily removed. 

The motion of a given moving element may be described 
relative to an infinite number of reference systems. However, 
these systems are not in general of the same significance. For 
let A, B, and C be three systems from which the motion of the 
moving element P may be observed. Suppose it is found that 
the motion of P relative to A is conditioned by that of B , but 
is independent of that of C. In such a case the motion of P is 
said to be related to />, which is known as a related reference 
system . (7, on the other hand, is an unrelated or ideal refer - 
ence system. Thus for the motion of a shot, the gun from which 
it is fired constitutes a related reference system. The velocity of 
a sound wave is determined, not by the motion of the source, 
but by the characteristics of the medium through which it passes. 
Hence in this case the source is an ideal reference body, while 
the medium is a related one. 

3. Principle of relativity. In the case of light, it has been 
generally recognized, ever since the vindication of the wave 
theory by Young and Fresnel, that the source does not consti- 
tute a related reference system. Recent analysis of the observed 
motion of certain double stars has confirmed this supposition. 
But most physicists have felt it necessary to postulate the exist- 
ence of an all-pervading medium in order to form a mental 
picture of the propagation of light waves through otherwise 
empty space. For a long time they were accustomed to attribute 
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to this medium, known as the ether, the properties of a related 
reference system. Finally Michelson devised an experimental 
method of measuring the velocity of the earth relative to the 
ether, based on the assumption that the ether is a related refer- 
ence system for the motion of light. Much to everyone’s 
surprise, this velocity turned out to be zero. Excluding the 
possibility of the earth’s being at rest relative to the ether, and 
one or two other equally improbable explanations, the only con- 
clusion to be reached was that the assumption that the ether is a 
related reference system for the motion of light was unjustified. 
The inference to be drawn from the result of this experiment, 
then, may be embodied in the following " principle of negation.” 

For the motion of an effect which travels through empty space , 
such as a light wave or one of the moving elements which form an 
electromagnetic or a gravitational field , there is no related refer - 
ence system . 

An immediate consequence is contained in the following 
statement. 

If a law governing physical phenomena which are conditioned 
solely by those effects which travel through empty space , is deter- 
mined from observations made in two different reference systems, 
the form of this law and the values of the constants entering into 
it can differ in the two cases only in so far as the geometry and 
devices for measuring time and distance, together with the units 
of these quantities, may differ in the two systems . Their relative 
motion in itself can affect neither the form of the law nor the 
values of the constants involved. This is the principle of general 
relativity. 

Consider two reference systems which have the same geometry, 
devices of the same character for measuring time and distance, 
and interchangeable units of these quantities. Such systems 
may be said to be reciprocal. It follows that 

A law governing physical phenomena which are conditioned 
solely by those effects which travel through empty space, has the 
same form and its constants have the same values for two mutually 
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reciprocal systems . In the subsequent discussion the phrase 
" principle of relativity ” will be understood to refer to this 
restricted form of the general principle. 

4. Reciprocal systems. Consider two reciprocal Euclidean 
systems S and S', such that all points of S' have the same 
constant velocity v relative to S. Let light travel in straight 
lines in S with a constant speed c. Then the principle of 
relativity requires that light shall travel in straight lines in S' 
with the same constant speed c . Let A and B be two points 
of either system a distance A r apart. Since the speed of light 
is the same in all directions, the time At taken by a light wave 
in passing from A to B, as measured in the system in which 
these two points are located, is the same as that taken by a 
light wave in travelling from B to A . Moreover, 

A r = cAt. (1) 

5. Differential transformations. Let a set of right-handed 
axes XYZ be fixed in S so that the X axis has the direction of 
the velocity of S'. Let a similar set of axes X'Y'Z', parallel 
to XYZ respectively, be fixed in S'. Let x , y , 2, and t be the 
coordinates of a point and the time at the point as measured 
in S, and x', y ', 2', and t' the corresponding quantities as meas- 
ured in S'. It is desired to obtain x', y', 2', and t' as functions 
of x , y, 2, and t. Let A' and B' be two neighboring points of S'. 
A light wave leaving A' at the time t arrives at B' at the time 
t 4 - dt x , and one leaving B' at the same time t reaches A' at the 
time t -f- dt 2 , these times being measured in S. If the coordi- 
nates of A' and B' relative to S at the time t are denoted by 
x, y, 2 and x 4 - dx, y + dy, 2 4 - dz respectively, the time dt! taken 
by the first wave to travel from A' to B' as measured in S' is 

dt ’ = fA dx + vdt ^ + % dy + Tz dz+ ft dt " 
and the equal time taken by the second wave in passing fro 
B' to A ' is 

. dt' , , N dt' . dt' , dt' 

dt' = — (- dx 4- vdQ - — dy — — dz + - - dt . 
dx cy dz ct 1 
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Subtracting, 

d ±{dx + \v (dt , + + + <* i - *«> = °- ( 2 ) 

and adding, 

<2/ = <?d£' = C L fa; 

But the squares of the distances travelled by the wave going 
from A! to B' and by that passing from B' to A' as measured 
in 8 are respectively 

<?dt\ = (dx + vdtj 1 + dy 1 + dz *, 

<?dtl = (dx - vdt^y + dy 2 + dz 2 . 


Hence 


<ftj— dt t — 


2v , 

(4) 

2 Vc 2 — D 2 ^ 

*2 .,2 rtr ’ 

(5) 


where dr‘=dx‘+ dy 2 + dz\ and v, is the component of v at right 
angles to dr. Substituting this value of dt,— dt 2 in (3), 

_i _^ + S!*V + ^ + 2:*=.o, 

1 — /3 2 \dx c dt / dv #3 




where /3 = -* 
c 

Now efc, d#, and dz are arbitrary. Hence their coefficients 
must vanish, that is, 

( 6 ) 

dx c dt 

fa 1 /"7\ 


^ = 0. 

dz 


Therefore the complete differential dt' is given by 
,, dt' u dt' dt' dt' . 
dt =H it + Tx* X + T/ y+ te' i ‘ 


Jl(it-Hdr). 
dt\ c / 
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To obtain dr', substitute (5) and (6) in (3). Then 

Giving dr the values dx v dy , and dz, it is seen that 

dx' = d ^dx v ( 10 ) 

dy> = 8 lVT=rFdy, ( 11 ) 

dz' JlLy/T^Wdz. (12) 

ot 

Now dx j in (10) is the distance of B' from A’ when the time 
is the same at the two points. If dx is the distance of the position 
of B' at the time t + dt from that of A 1 at the time t, 


So (10) becomes 


dx x = dx — vdt. 

dt' 

dx' — — (dx — vdt). 


Since the units of length in S and S' are interchangeable, it 
follows from symmetry that dy' = dy, and dz' = dz. Hence if 

k = l= , 

V1-/8 2 


Hence the differential transformations between the two sys- 

terns are / a \ / R \ 

dt 1 = k\dt dx\ or dt — k ydt f 4- — dx'J, (13) 

dx' = k (dx — vdt), dx = k (dx' + vdt'), (14) 

dy' = dy, dy — dy', (15) 

dz' = dz, dz = dz) (16) 
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where the second column is obtained from the first by solving 
for the unprimed differentials, or by changing the sign of v. 
From these expressions it follows that 

dr ' 2 — c 2 dt ,2 = dr 2 — c 2 dt 2 

is an invariant of the transformation. 

Now the velocity v f of a point of S relative to S r is obtained 
by dividing dx f by dt\ x remaining constant. Therefore 

v'— — Vy 

as might have been expected from considerations of symmetry. 

6. Space and time transformations. Integrating the differ- 
ential relations (13) to (16), and determining the constants of 
integration on the assumption that the origins of the two sys- 
tems are in coincidence when the time at each is zero, 


SQJ 

1 

Jl 

so 

or t — k(t r + — x'\ 

(17) 

\ e / 

\ c 1 

x! — k(x — Vt)y 

x = lc(x! - j- Vt r> )y 

(18) 

y'=y> 

y=y'i 

(19) 

z'= Zy 

z = z r . 

(20) 


These four relations are known as the Lorentz-Einstein 
transformations. 

Consider a moving element whose velocity components rela- 
tive to S are V x , V y , and V z . From the differential relations of 
the preceding section it follows at once that 



( 28 ) 
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For the resultant velocity 


F' 2 = c 2 



Hence if V= e, V'= c also, as should be. Now 



so the velocity of light is the only velocity independent of v 
which is the same for the two systems. 

Suppose that the velocity of S r relative to S is nearly as 
great as the velocity of light. Then f3 = 1 — 8, where 8 is small. 
Consider a body moving in the X ' direction with a velocity 
relative to only slightly less than the velocity of light. 
Then F r '=c*(l — e), where e is small. Equation (21) gives for 
the velocity of this body relative to S, 


V = c 


2 — 8 — 1 


2 — 8 — e + Se 



whence V x is less than c. For example, if v = 0.9 c , and Fj = 0.9 c , 
V. x would be 1.8 c according to nineteenth-century conceptions 
of space and time. But the addition theorem of velocities just 
obtained from the principle of relativity gives F x = 0.994 c. 

The relations between components of acceleration as measured 
on the two systems are obtained by differentiating equatior|s (21) 
to (28) \yith respect to the time, remembering that 

d dt d 
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f.= 


These relations are 

fx 


F 


MW 


or 


/*=• 


./ as 


F 


MW 


( 24 ) 


fi 




f y +-(fx-m 

/= f (25) 

•'* / t//\3 


f = ' 

J V 


MW 


F 


MW 


/; = 


c 


F 


mw 


1 — — — < 26 ) 


F( 1+ *sy 


Suppose that the moving element is at rest in Then the 
second column becomes 

■iV 

(27) 


f 

J? */ as 


f = ‘ / -K. 
•' y 


fi, 

F 




F* 


(28) 

(29) 


7. Four-dimensional representation. The Lorentz-Einstein 
transformations can be represented very simply by a rotation in 
a four-dimensional manifold. For consider a set of rectangular 
axes XYZL in four-dimensional space, such that the distances 
x , y, z of a moving element from the origin of S are measured 
along the first three axes, and the quantity l = ict along the 
fourth axis, where i = V— 1. The position of a moving element 
at a given time is represented by a point in this space, and the 
locus of the positions of such an element at successive instants 
by a line. This line is called the world line of the moving 
element. Thus the world line of a body permanently at rest 
relative to $ is a straight line parallel to the L axis. The world 
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line of a point moving with velocity V relative to S is inclined 
to the L axis by an angle <f> such that 

tan (f> = — i — » 

c 

Hence the world line of a point of S' is parallel to the XL 
plane and makes with the L axis an angle a given by 

tan a = — i/3, sin a — — ik/3, cos a = k. 

Therefore in terms of x , y, z, and l the Lorentz-Einstein trans- 
formations take the form 


• V = l cos a + x sin a , 

or l — V cos a — x f sin a. 

(30) 

x r =x cos a — l sin a , 

x = x' cos a + V sin a , 

(31) 


y = y\ 

(32) 

z r = z. 

Z — z r . 

(33) 

amounting formally to a rotation of the X and L axes 

through 


the imaginary angle a . 



CHAPTER II 


THE RETARDED FIELD OF A POINT CHARGE 


8. Electric field. Continuous lines may be imagined to spread 
out from every elementary electric charge in such a way as to 
diverge uniformly in all directions when viewed from the system 
in which, at the instant considered, the charge is at rest. These 
lines are called lines of force and, taken 
together, they constitute the charge’s field. 

The number of lines emanating from an 
element of charge de will be supposed to 
be very large, no matter how small de 
may be. A bundle of M lines, where M 
is a very large number, will be considered 
to constitute a tube of force. The field 
strength, or electric intensity , E, at a point 
in a field, is a vector having the direction 
of the lines of force at that point and equal in magnitude to 
the number of tubes per unit cross section. Thus if dN tubes 
pass through a small surface of area da whose normal makes 
an angle 6 (Fig. 1) with the field, the magnitude of the field 
strength is given by 

JE — • 

da cos 6 



Hence the component of the field strength parallel to da is 




dN 

da' 


or, in vector notation, 


r> dN 1 

*»= i?*"- 


The ratio of the number of tubes of force diverging from a 
charge to the magnitude of the charge determines the unit of 

20 
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charge. The simplest unit, and the one which will be used in 
the following pages, is that which makes the charge at any point 
equal to the number of tubes of force diverging from that point. 
This unit is the one advocated by Heaviside and Lorentz and, 
as will appear later, is smaller than the usual electrostatic unit 
by the factor ^ 

V4 7 r 

9. Motion of a field. Consider an electric field which is being 
observed from the two reference systems S and S r of the pre- 
vious chapter. The principle of relativity requires that the 
velocity of the moving elements comprising the field shall have 
the same numerical value no matter whether observations are 
carried on in S , S', or some other system reciprocal to one of 
these. In section 6 it was shown that the velocity of light is the 
only velocity which satisfies this condition. Hence the moving 
elements constituting an electric field must have the velocity 
of light. 

Suppose a charged particle to be permanently at rest in S. 
Although the moving elements constituting its field are in 
motion with the velocity of light, the lines of force themselves 
are stationary. Hence the motion must be entirely along these 
lines. Now consider a charged particle moving with a constant 
velocity V relative to S. As the charge carries its field along 
with it, the velocity of a moving element will be along the lines 
of force only at points in the line of motion. At all other points 
this velocity will make an angle with the lines of force which will 
be greater the greater the speed. Therefore in general the com- 
plete specification of an electric field due to a charged particle 
requires the knowledge at every point of the values of two vec- 
tors, the field strength E and the velocity c. Both magnitude 
and direction of E must be given, but as the magnitude of c 
is known its direction only is required. 

10. Transformation equations. Suppose that the field strength 
E and velocity c relative to S are known for the field due to a 
point charge. It is desired to find the values of these quantities 
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as measured in S'. Let P and Q (Fig. 2) be two neighboring 
points on the same line of force, the coordinates of Q rela- 
tive to P being dx, dy , dz when the time is the same at 
the two points 
according to 
the standards 
of S. Then if 
the time at Q 
should be later 
by dt than that 
at P, the coor- 
dinates of Q 
relative to P 
would become 

dx x — dx -b c x dt, 
dy=dy+c y dt, 
dz y — dz + c z dt . Fig. 2 

Now, in order that the times at P and Q 1 should be the 
same when measured in S r , equation (13) of section 5 shows 
it to be necessary that 

dt = —dx,. 
c 

Hence when dt! — 0, 

, dx 

dx. = 

1 — /3— 

o 



(3 C -±dx 

d ?J l =dy + 7» 

1— /3 — 


dz L = dz -}■ 


/3^dz 

c 

1 — /3 — 
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Substituting in the transformation equations (14) to (16), 
second column, of section 5, 

= , , 


:(l — /3 S 


dy f = dy + * 


/3 ^ dx 
o 


dz'= dz + 


1-/3^ 

c 

p- ±dx 

c 


are found to be the coordinates of Q relative to P as measured 
in S' when the time is the same at the two points according to 
the standards of this system. The position of Q at this time is 
shown by Q ! on the figure. Consequently the lines of force in 
S r extend from P to Q\ instead of from P to Q as in & 

Let hN be the number of tubes of force passing through the 
area SySz in S. Then, since Sy f = Sy, and Sz' = Sz , 

bn sir _ 

x 8/Sz'. SySz 


Moreover, 


Hence 


and similarly, 


dx dy dz 

dx f dy' dz! 

dy — /? (~dy — ^ dx'j 


El = JcE a 


= |(ex E)t }’ 

E’ t = k + ^(cx E) v } . 
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The magnetic intensity H is defined by 


H = -(cxE). 

c 



Therefore E' x — E x , 

r "1 ^ 

a) 

E’ u =k{E y -/3lf z } = k 

(^+-( vxH X)’ 

(2) 

E< = k{E z + pH y }=k\ 

+ 

! M> 

X 

w 

Jr' 

(3) 


where the primed and unprimed quantities may be interchanged, 
provided the sign of v is changed. 

The transformations for c are obtained at once from the equa- 
tions (21), (22), and (23) of section 6 for transforming velocity 
components. They are 



From the transformations for the components of E and c 
those for the components of H are readily deduced. For 


c' r f 

//'= E r — E r 
x c z c y 


by definition. Substituting the values of the components of E' 
and c f in terms of the unprimed quantities in this identity and 
in the corresponding expressions for H y and H z , it is found that 


H' v = k{H v +!3E z } = k{ 

r 

fci 

1 

ft 1 ^ 

< 

X 

3 

ts 

/-N 
l>- 00 

jn=k{ir-pE v }=k. 

k-|(vxE) z }' 

(9) 


A field formed by the superposition of the individual fields 
of a number of charged particles is termed a complex field, in 
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contradistinction to the simple field of a single elementary charge, 
and the electric and magnetic intensities in such a field are defined 
as the vector sums of the corresponding intensities of the com- 
ponent simple fields. It is to be noted, however, that the velocity 
c, in so far as its direction is concerned, always refers to the field 
of an element of charge, and never to the resultant of a num- 
ber of such simple fields superposed. If it is desired to avoid 
explicit reference to the components of a complex field, the field 
must be described by means of equations which do not involve 
the direction of motion of the constituent moving elements. 

Since the transformations that have been obtained for E and 
H due to a single point charge are linear in these quantities, 
they apply as well to complex as to simple fields. 

11. Point charge at rest. Let e (Fig. 3) be a point charge 
momentarily at rest at the origin of S at the time 0. It is 
desired to find the 
field strength E at 
P at a time t = r/c, 
where r is the dis- 
tance of P from 0 . 

Two moving elements 
leaving e in slightly 
different directions at 

the time 0 will be at fa Flo 3 

P and Q at the time t. 

At a time dt, e will still be at 0 (to the second order of small 
quantities). Consequently a moving element coming from e at 
this time and belonging to the same line of force as that at P 
will reach some such point as A by the time t. Similarly, one 
belonging to the same line of force as that at Q will reach B by 
the time t. Hence, if r x denotes the distance OA, 

r = ct, 

and r x = (c + dc ) (t — df) 

= ct — c dt -f- td c, 

where |c + <Zc| = |c| necessarily, or dc is perpendicular to c. 
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Put p = r - r r 

Then p = (cdt — ^ dcj , 

p.r = crdt . 

Now since the lines of force, as viewed from S, diverge uni- 
formly in all directions from 0 at the time 0, the number of 
tubes per unit area passing through a small surface at P with 
normal parallel to r is 

e 

4 7 rr 2 


Hence 


4 7 rr 2 p cos QPH 

f P 

4 7i rr p-r 

e f 

4 irr 2 c \ c dtJ 


( 10 ) 


Suppose that the charge e has an acceleration f relative to S. 
Then at the time dt it will be at rest in some reciprocal sys- 
tem which has a velocity f dt relative to S. As the lines of 
force diverge uniformly from the charge when viewed from the 
system in which, for the instant, the charge is at rest, and as 
the velocities of the moving elements constituting the portions 
of these lines in the immediate vicinity of the charge are along 
the lines themselves when observed from this system, it follows 
that if two moving elements, one of which leaves e at the time 0 
and the other at the time dt , are to lie on the same line of force, 
the velocity of the second must make the same angle in S f with 
the direction of f as that of the first does in S. If the velocities 
of these two elements are denoted by c f and c, and if the X f and X 
axes are taken parallel to f, 




= c 


( c x + d cj) — fdt 

c x fdt 


But, from (4), 


1 
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or 


Let a (Fig. 4) be the angle between c and f. Then 
c x = c cos a, 
dc x = — c sin ada. 


Hence 


Now 



which becomes, in vector notation, 

<? 8 


dt. 


(ii) 


Substituting in (10) it is seen that if a charge e which has an 
acceleration f is momentarily at rest at 0 at a time 0, the field 
strength at a point P distant r from 0 at a time r/c is given by 

*'>*']]• (12) 

where the heavy brackets are used to denote the fact that the 
quantities contained therein are retarded ; that is, these quan- 
tities refer to the effective position of the charge, or its position 
at a time r/c earlier than that for which the field strength at 
P is to be determined. 

Since H = -(c x E) 

c 

by definition, H £{(* x c) x c} x cjj . (IS) 

12. Point charge in motion. Consider a point charge e which 
is passing the origin of S at the time 0 with velocity v and 
acceleration f. Choose axes so that v is along the X axis. Then 
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this charge is at rest, at the instant considered, in a System S' 
which has a velocity v in the X direction relative to S. So 


and 


E, = 4^[ C ' + ? (f><C, H 

*c[“ x< ^ xcl * xc> ] 


H' = 


4 7 rr r2 


are the values of E f and H f at a point P distant / from 0 at the 
time r r /c . It is desired to determine E and H at P at this same 
instant. Since the velocity of light is the same in the two systems, 
the time at P will be r/c in S when it is r r /c in S f , Hence the 
result of the transformation about to be carried through will 
give E and H at P at the time r/c. 

Let a ' be the angle in S r which the line OP makes with the X r 
axis. Without loss of generality this line may be supposed to lie 
in theX'F' plane. The Lorentz-Einstein transformations give 


whence 


r' cos a ' = hr (cos a — /3), 


r' sin a f = r sin a ; 


r l = kr (1 — /3 cos #), 


, cos a — /3 
cos # = : 


1 — /3 cos a 


sin a 1 


sin# 


Jc,(l — /3cos #) 


From (4), (5), and (6> c' x = , 

o' = S' , 

r k(l — /3 cos a) 

4 = 0 ; 

and from (27), (28), and (29) of section 6, 

fl = W v , 
fl = W, 
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Now 
E, — EL 


= 4^J' i+ £ K '' X ‘' )XC '4 

E=Tc{El+m 

x c ’> x c ’>'{ 1+ ^7}-« f x ^ xc ’>.^)]' 

E t = k{E' z -pH> y } 

=45v[ 0 + ?« f ' xc ’) x ^{ 1 + ^7}]- 

Substituting for the primed letters their values in terms of 
the unprimed ones, and reducing, 

^ r vie 2 ~1 

E * = 4 TrrVc (1-/3 cos ay[ Cx ~ V + ~7 X X C ^J ’ 

e r t 

Ey= 47rr 2 Fc(l-/3cosa) 8 L C ' ,4 ”^^ f X (C- ^ } X c \|’ 

^- 4^a-^o Sg y[° + ?« f x (‘-'M x 4 


Jg _ 

* 4 irr 2 k 2 c (1 — ft cos 


4 irr 2 k 2 c (1 — /3 cos a ) 


4 7rr 2 Jc 2 e (1 — /3 cos a ) 


4 7 rr 2 k 2 c 


i r 


c — v H -g- {f x (c v)} x c I. (14) 


H may be obtained from H r and E r in an exactly analogous 
manner, but it is more easily found directly from the above 
expression for E. For 1 

He-(cxE), 

whence c 


e 

4 7rr 2 kV^l- 


s[~ cx 


v +i c X({ f *(«-▼)} 

0 


xc)] . 


These expressions for E and H, it must be remembered, give 
the values of these vectors at P at the time rjc in terms of v 
and f at 0 at the time 0, c having the direction of the line OP. 
In other words, all the quantities within the heavy brackets are 
retarded. Each expression consists of a part involving the inverse 
second power of the radius vector r, and a part involving the 
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inverse first power only. The latter depends upon the acceler- 
ation of the element of charge, and the part of the field which 
it determines is known as the charge’s radiation field, 

13. Retarded potentials. In differentiating expressions such 
as those involved in (14) and (15) account must be taken of the 
fact that the quantities enclosed in the heavy brackets are 
retarded. Let [ fir ] be a scalar whose value at P at a time t is 
given in terms of the position and velocity of a charged particle 
at 0 at a time t — r/c, where r is the distance OP. Then, if the 
coordinates x , y, z of P relative to O remain unchanged, 

■w- BH- 

r 

where, since dt = [dt] + ~ — - 


it follows that 


Consequently, 


or, symbolically, 



(16) 


On the other hand, if x changes by an amount dx while y, z, 
and t remain constant, 

where [dx] = dx , 



m=- 
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Hence 


dx 


K] 


c.y 




and similar expressions hold for differentiation with respect to 
y and z. Therefore 


V = [V] — 


1- 


c.v 


(17) 


Consider the retarded scalar potential function 

M 




It will be shown that E and H due to a point charge may be 
expressed as derivatives of this function in much the same way 
as the field strength due to a point charge in electrostatics may be 
given in terms of the gradient of the potential function 1/r. For 



vm= 


and 




'-7K 

c-f| 



r 2 c(l 

c-v\ 8 




e 2 / 




c{-(l-/3 2 )- 

7 c -'} + T ( l 

c*v \” 

C 2 ) 


r 2 o(^ 

i_±iY 






J 


c a 

1 

{(l-/3 2 )-( 

‘-t) 

+ 5 c - f }y 


r 2 c(^ 

1-- 

.v^ 

T) 

8 

1 


Therefore 

1 d 




(l-/3 2 )(c-v) + Jjfx(c-v)j xc 
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and comparison with (14) shows that 
Moreover, 


E = — V hH — — 

4 7T VlTJ 4 7 TO St 


M 


Vx[lf].[f]Vxp]-[']xVW 

’ r 



and comparison with (15) shows that 


H = f-V 

4 7 r 




Whether the field is due to a single point charge or to a 
number of such charges, if 


-2 


w 


e 


(18) 


and 


-2 


e 

4 7 r 


D»] 


then 


(19) 

(20) 
( 21 ) 

Hence an electromagnetic field may be specified by the values 
at all points and times of either the two vectors E and H or 
the scalar potential <j> and the vector potential a. 


E = — V<£ — ~a, 
c 

H = V x a. 



CHAPTER III 


THE SIMULTANEOUS FIELD OF A POINT CHARGE 

14. Constant velocity. Let a point charge e which has a 
constant velocity v relative to S be at the origin O of S at the 
time 0. It is desired to find the values of E and H at a point P 
(Fig. 5) at the same instant in terms of 
the coordinates of P and the velocity v. 

Choose axes so that the velocity of the 
charge is along the X axis and the point 
P is in the X Y plane. The point Q occu- 

v 

pied by the charge at the time is its 

effective position. Hence equation (14), 

section 12, gives for the field strength q g 1 — 

at P at the time 0, Fig. 5 

e l (c — v} 

4 7 rr 2 k 2 c (1 — /3 cos a:) 8 

Now QO=fir,* 

and therefore the vector c — v, and consequently E, have the 
direction OP. 

From the geometry of the figure 

r 2 = /3V 2 + fi? 4- 2 /3rP cos 0, 

or — = Vl— /9 2 sin 2 0 — /3 cos0. 

r 

Hence -|c — v| = Vl— /3 2 sin 2 #— j3 cos#. 

c 

Moreover, 

1 — /Scos a = Vl — /3 2 sin 2 # {V 1 — /3 2 sin 2 # — /3 cos #} • 

* In order to avoid cumbersbme equations in this and the succeeding article, 
the retarded quantity r is not enclosed in brackets. 

33 
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Therefore E = — — — — — , • (1) 

4 ttH* (1— 


The magnetic intensity H is directed upward at right angles 
to the plane of the figure. Its magnitude is given by 
H=E sin(0-a). 

But csin(0— a)=vsmO. 

e (1 — /3 2 )/3sinfl 


Hence 




4 irli 2 (1-/3W0) 1 


( 2 ) 


The expression for E shows that the lines of force diverge 
radially from the moving charge (Fig. 6), but, instead of spread- 
ing out uniformly in all directions, as in the case of a static 
charge, they are crowded together in 
the equatorial belt and spread apart 
in the polar regions. The greater 
the speed the more pronounced this 
•disparity, until, if the velocity of 
light is attained, the entire field is 
confined to the equatorial plane. 

If lines are drawn so as to have 
everywhere the direction of the vector 
H, these magnetic lines of force will 
be circles in planes at right angles 
to the line of motion with centers lying on this line. If the 
magnetic lines, like the electric lines, indicate by their density 
the magnitude of the corresponding vector, a similar crowding 
together in the equatorial belt and spreading apart in the polar 
regions will exist. However, the total number of magnetic lines 
of force in the field, unlike that of the electric type, will be 
greater the greater the speed of the charged particle. 

15. Constant acceleration. Consider a point charge moving 
with an acceleration <j>, which always has the same value relative 
to that system, reciprocal to S ', in which the charge happens to . 
be at rest at the instant considered. Let this charge come to 
rest momentarily at the origin 0 of S at the time 0. It is desired 



Fig. 6 
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to find the values of E and H at a point P (Fig. 7) at the same 
instant in terms of the coordinates of P and the acceleration <|>. 
Choose axes so that <(> is along the X axis and P lies in the 
XY plane. The effective position of the 
charge is the point Q y which it occupied 

at the time — • From (27), section 6, 
c 

/=<K 1-0*) 4 - 

Integrating, the velocity at Q is found 
to be given by 


0 = 


<j>r 

1? 




1 + 


<*>V 



and the distance OQ by 
OQ- 




If the components of E along and at right angles to QP are 
determined separately, it follows from equation (14), section 12, 
that 

e 1 


P<t — o, 


[ 


4 7T7*% 2 (1 + /3 cos a) 




(3) 


showing that E is parallel to c. 
Now, from geometry, 


= E - 


N 


l+^rCO80 + 


4 c 4 


1 + ^T COS0 

c 2 


jBcos0=rcosa + OQ. 

1 + ^t cos# 

c 




j>y 

c 4 


Hence l+j3cosa 



36 AN INTRODUCTION TO ELECTRODYNAMICS 


Substituting in the expression for E r . 


1 + vj± cose + ' 


Er 4 ttE 1 


To find the components of E along and perpendicular to OP 
it is necessary to obtain the values of cos 8 and sin 8. From the 
geometry of the figure it follows that 


cos 8 = . 


1 + cos 6 
2 if 

\ , <t>R a 7 

1 + cos 6 + ~—t 


whence 




Er 4 wJF 


cos 6 -f - 


1 + s 7 cm# 


1 + ^COS^ 




i* c™«+^ 


To obtain the equation of the lines of force, use may be made 

of the relation ^0 

E ■— = tan 8 
dE 

<j>E . a 


l + ||cos^ 

The solution of this differential equation is 

cotd + P^-X -~b, 

2 cr sm 6 
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where b is the constant of integration, 
the form 


Writing this equation in 


(^R cos 6 + + QR sin 0 — by 2 


: b Z + 




( 7 ) 


it is seen that the lines of force are circles passing through 0 
with centers in a plane at right angles to the X axis and at 
a distance —c 2 /<f> from 0. Fig. 8 shows the section of the field 
cut by a plane through the* 

X axis. The full lines repre- 
sent the lines of force. 

Consider an electric field 
all parts of which happen to 
be at rest in a single system 
S at the same time. A sur- 
face in $, so constructed as to 
be everywhere normal to the 
lines of force, is called a level 
surface . For a field perma- 
nently at rest in S level surfaces are identical with the equipo- 
tential surfaces of electrostatics. In the case under discussion of 
a charged particle moving with constant acceleration, E has been 
shown to be parallel to c simultaneously at all points. Hence all 
parts of the field are at rest at the same time and consequently 
level surfaces can be constructed. The differential equation of 



Fig. 8 


these surfaces is, 


J? dO * 

R — = — cot o 
dR 




of which the solution is 


4>R 

2 > 


sin# 


- cos0— — + sin0} 2 == (8) 

where h is the constant of integration. This is the equation of 
a family of spheres with centers at the effective positions of the 
moving charge. Their traces are shown by broken lines in Fig. 8. 

Since E is parallel to c at the instant considered, H is every- 
where zero. 
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16. General case. The retarded expressions for E and H 
deduced in the preceding chapter show that the field at a point 

P and time 0 is conditioned by the motion of the charge pro- 

M 

ducing that field at a time — ^ > where [r] is the distance of P 

from the effective position of the moving particle. Therefore 
the specification of the entire field at the time 0 involves the 
complete past history of the charged particle, Since, for physical 
reasons, the motion of this particle must be continuous, the past 
history of its motion is contained in its present position, veloc- 
ity, acceleration, and higher time derivatives of the positional 
vector. Hence the simultaneous values of E and H may be 
expressed as series in these quantities. While these series may 
fail to converge for distant portions of the field, or for very 
rapidly changing motion, their form will make evident their very 
rapid convergence for all cases to which they will be applied. 

If the point at which E and H are to be evaluated is chosen 
as origin, and if [x], [y], [ z ] are the coordinates of the effec- 
tive position of the charged particle relative to this origin, 

“--DO- 

etc., and the x component of the retarded expression (14), 
section 12, for the field strength, may be written 

If Xy y , z are the coordinates of the simultaneous position of the 
charge, and v, f, f • • • its simultaneous velocity, acceleration, etc., 
then ri -| r 21 i i .. r*. 4 i i ...f». 6 i 


\x\-x V x c +2 Jx f qJx c a +2 ¥ x c* 


j -f-n 

12 Q Jx c 6 


r i f M , 1 i M 

[*y = v x —Jx ~ + o fx^-T- 


V x c 8 ' r 24* / * c* ' 


lJx\ Jx Jx q I 2 gZ *’ 

and similar expressions hold for the y and z components. Put 

m -X a - V X „ -f* r S _/y , - f/ y = fx r * 

m x = -, /3 x = -, yx =—, S x =—, ex= x—, 
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Then ^ s 

[r 2 ] = r 2 {l- 2 p.m + (ym + /3 2 ) & J - g (8-m + 3 vP) ^ 
+ i( € .m+48.p + 3 7 2 )M 
-i(^m+5 € .p+108.v)M...}. 

[r] 

Now put T = — , 


a=p.m&, 

Jsym?, 

6* == (8*m -f B 

d == (e*m + 4 8*p + 3 v 2 )& 4 , 

0 ee (£.m + 5 €*p + 10 8*-y) £ 5 , 
whence the previous equation becomes simply 

t 2 = 1— 2ar4-&T 2 — ^cr 8 -*-^ dT* — ^eT b * (10) 

aild t 1 + = + 0 + 1 ct8 - n dr4+ h eTb " '} 

= j 
_ Fr 2 ‘ 


After some reduction it is found that 





L{,»,(i- 2 »t + 0 +|«t- 


5 , , 6 , 

i2* + 40" 



+ o 

-I^(-i + 0 +0+ l CT ' _ n* A+ ^ OT ****) 

■ + 1 s ^ v (- 2 + aT + 0 + \ - ^ T< + ^ ct6 • • •) 

3 + 2 ar + 0 + 0 - ^ *- 4 +^er*. . ) 

+ M^ T6 (~ 4 + 3aT + °~^ S+0 + ^ eTB ''') , ‘^ 

: Ft 8 ’ 
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Hence -^—krI~\L (11) 

* 4 irr 2 v J 

Returning to (10) and solving for r by successive approxi- 
mations, it is found that 

t=1— a(\ — | a -h 0 -f ^ + 0) 

+ |-6(l-2a + fa 2 +0) + |5 2 (l-f a) 

— S c (l“^ a ~f-4 # 2 ) — j hc+ ^ d (1— 4 a) — j%Q e •• •. (12) 


Substituting this value of t in (11) and reducing, the fol- 
lowing expression is obtained for the X component of the 
electric intensity, 



The X component of the magnetic intensity may be obtained 
most easily from the relation 



which gives, after some reduction, 

H * = ~ 4 ^ ( X “ 1 “ 2+ \ h • ' •) 

+ ak(y y rn z - 7 a wi„) “ \ k XK m *~ K m v)+ ^ ^(e,/n s - e a m y ) 

-|^(7A-yA)+|**C*A- W • • < 14 ) 

Equations (13) and (14) give E and H at the origin in terms 
of the simultaneous coordinates, velocity, acceleration, etc. of 
the charge e. To obtain the values of E and H at a point #, y, z 
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due to a charge e at the origin, the signs of the coordinates in 
these two equations must be changed, giving 


JBL 


elc 


4 irr 1 




3 2 , 15 4 

2 a ' + T a ! 


1, ^9 „ 

2^1 + 4 a i 


■ + S*f 


+ i^ 3 ( 1 +|« I -2/ V . •)-|^ 4 (l+!« 1 . . .)+± . • •}, (15) 


- ajc (j1y m z - r iz m y) ih m z - S * m v) ( e « m z - «*«*„) 

+ A- 7A>- f*W*- W ■ • •}> (16) 

where = p-m/c, 

?> L = y*mk 2 , 

c \ ( 6 *m — 3 -y * P) /c 8 , 

^ = ( € .m - 4 8-p - 3 7 2 ) & 4 , 
e 1 == (£*m — 5 €*p — 1 0 8-7) & 5 . 



CHAPTER IV 


THE DYNAMICAL EQUATION OF AN ELECTRON 

17. Electrical theory of matter. All matter will be assumed 
to be made up of positive and negative electrons. An electron 
will be defined as an invariable charge, of magnitude approxi- 

matel y 4.77(10)- lo V4^ 

Heaviside-Lorentz units, distributed over a surface which is 
spherical in form to an observer in that system, reciprocal to S, 
in which the electron happens to be momentarily at rest. A 
positive electron will be considered to differ from a negative one 
only in the sign of the charge involved and the radius of the 
spherical surface over which it is distributed. 

The electromagnetic force dK on an element of charge de, as 
measured in that system, reciprocal to S, in which this charge 
happens to be at rest at the instant considered, is defined as the 
product of the field strength E by the charge de. The extension 
of this definition to the case of a system in which de is not at 
rest, will be given in the next section. 

The distribution of charge on the surface of an electron will 
be supposed to be such as to make the tangential force due to 
its own field zero at all points of the surface, all measurements 
being made in that system relative to which the electron is 
momentarily at rest. This assumption is introduced merely for 
the purpose of simplifying the analysis (section 20) involved in 
determining the dynamical equation of an electron moving with 
constant acceleration. To the number of terms to which the 
analysis is carried ip the general case (section 21), no change 
in the dynamical equation is introduced if this hypothesis is 
replaced by the more probable assumption that the distribution 
of charge is such as to make the tangential force due to the 

42 
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total field, that is, the resultant of the impressed field and the 
electron’s field, equal to zero, or by the simple assumption that 
the distribution of charge is always uniform. 

18. Dynamical assumption. In the previous chapters the 
discussion has been concerned with the determination of the 
field of a charged particle. It must be borne in mind, however, 
that the lines of force constituting such a field are nothing 
more than convenient geometrical representations to be em- 
ployed in describing the effect of one charged particle on 
another, and that no reason exists for attributing a greater 
substantiality to them than to any other arbitrary convention, 
such as, for instance, parallels of latitude on the earth’s sur- 
face. The representation of a field by lines of force has led to 
the concept of electric intensity, and the electromagnetic force 
on an element of charge, as measured in the system in which 
the charge is momentarily at rest, has been defined in terms 
of this quantity. In order to pass from these definitions to 
the quantitative description of the effect of one electron on 
another, it is necessary to introduce the following dynamical 
assumption : 

The motion of an electron is such as to make the total electro- 
magnetic force on it , as measured in that system, reciprocal to S, 
in which it happens to be momentarily at rest , equal to zero . 
By the total electromagnetic force is * to be understood the 
resultant of the force due to the impressed field and that due 
to the charge’s own field. With forces which are not electrical 
in nature, such as must exist if a dynamical explanation of the 
stability of the electron is possible, the present discussion is 
not concerned. While extra-electrical stresses on a single elec- 
tron may be of great intensity, their resultant will be assumed 
to be always zero. Moreover, such forces will he supposed to 
be comparatively negligible when the effect of one electron on 
another is under consideration. Thus no account will be taken 
of the gravitational attraction between two electrons, as it will be 
deemed quite unimportant compared to the electrical attraction 
or repulsion. 
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Consider an electron which is at rest in S r at the instant 
considered. Then the dynamical assumption stated above re- 
quires that r 

H / Wde = 0, (1) 


Substituting for the components of E r their values in terms 
of the components of E and H, as given in equations (1), (2), (3), 
section 10, it is found that 

^ EM = 0 , 




j"{ 7 ^+^(vxH) ;/ j <fc = 0, 
f{^+*(vxH),}<fc = 0. 


Hence it is natural to extend the definition of electromagnetic 
force given in the preceding section so as to read : 

The electromagnetic force d K on an element of charge de, as 
measured in a system, reciprocal to S, relative to which the 
charge has the velocity v at the instant considered, is defined by 

iK = |E + i(vxH)|c?e. (2) 


Then the dynamical assumption may be stated in the more 
general form: 

The motion of an electron is 'such as to make the total electro- 
magnetic force on it , us measured in any system reciprocal to S, 
equal to zero . Thus the dynamical equation of an electron may 
be found directly for any system, no matter whether the elec- 
tron is at rest in that system or in motion with respect to it. 
However, in order to avoid unnecessary analysis, the method 
pursued will be first to deduce this equation relative to that 
system in which the electron is momentarily at rest, and then 
to extend it to other systems by means of the transformations 
already derived. 

19. Constant velocity. Consider an electron permanently at 
rest dn S r . Relative to an observer in this system the electron is 
a uniformly charged spherical surface of radius a with a uniform 
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radial field. To an observer in S however, this electron has 
a constant velocity v along the X axis, and the transforma- 
tion equation (18), section 6, shows that its dimensions in the 
direction of motion are shorter in the ratio 1 : k when viewed 
from this system, while those at right angles to this direction 
are unchanged. Hence to an observer relative to whom an 
electron is moving its surface is that of an oblate spheroid with 
the short axis in the direction of motion. 

Describe two right circular cones with vertices at the center O' 
of the electron and axes along the X axis such that elements of 
the cones make angles 0' and 0 r +d0 r respectively with their 
common axis. If e is the charge on the electron, the number of 
tubes included between the cones is 


Hence 


COS0 

COS 0 = — 7~::~ — y 

Vl — /S 2 sin‘ 2 0 

. (l-/3 2 )sin0cZ0 

(1-0W0)* 

_ e (1 — /3 2 )sin 0d0 
dX=- ± r , 

2 (l-/9‘ 2 sin 2 0)* 


x ■ y 

and the electric intensity in S at any distance E from O' is 

(Ml 

4 ' Ji (l-/3’sin'«) 


The magnetic intensity is given by 

#=-|cxE| 

c . 

= Efi sin 0, 

id substituting for E its value from (3), 

sin l 

4 Trie 2 (l-^sin^) 1 


( 4 ) 
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Comparison of these expressions with (1) and (2), section 14, 
shows that the external field of an electron moving with con- 
stant velocity is the same as that of an equal charge located at 
its center. It follows from symmetry that the resultant force 
on the electron due to its own field is zero. Hence the dynamical 
assumption requires that the impressed force shall be zero as well. 

20. Constant acceleration. Consider an electron each point of 
which moves with an acceleration which always has the same 
value relative to that system, reciprocal to S , in which this point 
happens to be at rest at the instant considered. Let <|> be the 
value of this acceleration for the point 0 of the electron. Choose 
axes so that <|> is along the X axis. Then (27), section 6, gives 
for the acceleration / of this point relative to S at any time 

/= *(1-/8*)*. (5) 


Integrating, the velocity of 0 relative to S is found to be 
given by 


fi = - 

and the displacement by 
x - 






C 

( 6 ) 



( 7 ) 


Consider a neighboring point P of the electron such that OP 
is parallel to the X axis and equal to d\ when 0 is at rest in S, 
Then, since O and P have constant accelerations relative to the 
systems, reciprocal to $, in which they happen to be at rest at 
the instant considered, the principle of relativity requires that 
the length OP as measured by an observer in any system recip- 
rocal to S when O is at rest in that system shall be the same as 
the length OP as measured by an observer in S when 0 is at 
rest in S. Hence OP = d\y/T^ (8) 


is the distance OP as measured in S when O has the velocity @c. 

Now, when 0 is at rest in S , P may have a velocity dp in 
the X direction. But the principle of relativity requires that dfi 
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shall be the same relative to any other system, reciprocal to S, 
at the instant when O happens to be at rest in that system. 
Hence, adding to d\ the difference between the displacements 
of P arid 0 in the time £, 

w = d\ + ^{l-y/T^W)+p c -^ (9) 

9 9 

is found to be the distance OP as measured in S when 0 lias 
the velocity f3c. 

Equating coefficients of like powers of /3 in the identical 
expressions (8) and (9), it is seen that 


dfJL = 0 , 



The first of these equations shows that when one point of the 
electron is at rest in S, every other point is likewise at rest . Inte- 
grating the second, , 

*— %r- <1»> 


1 + 4 

c 


where <f> 0 is the acceleration of 0. This equation shows that points 
on the forward side of the electron have smaller accelerations 
than those on the rear. Such a difference is obviously necessary 
in order to produce the progressive contraction of the electron 
required by the principle of relativity as its velocity relative 
to 8 increases. 

Obviously, the relations just obtained between the velocities 
and accelerations of points of the electron under consider- 
ation apply equally well to points of the field of Fig. 8, p. 37. 
Hence any one of the level surfaces of this field, such as that 
upon which the point P lies, may be considered to constitute 
the surface of the electron. As the charge is distributed entirely 
on this surface, it is necessary, in order that the external field 
should be the same as that due to an equal charge at 0, that 
the density of charge should be everywhere equal to the electric 
intensity just outside this surface. 
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If e is the electron’s charge, and the strength of the external 
field, the impressed force is 


K i =«Ef 


( 11 ) 


In Fig. 8, Q is the geometrical center of the electron and a 
its radius. The electric intensity E 2 at its surface due to an 
equal charge at 0 is given by 




4 7 raV (1 -f ft cos af 


( 12 ) 


where 


fir 






(f> 0 being the acceleration of the point 0 . As the field due to the 
electron vanishes everywhere within its interior, the resultant 
force K 2 on this charged particle due to its own field is 

K 2 = t r E 2 cos ada, 


where dcr is an element of the surface. Substituting for E its 
value in terms of a and integrating, 


k 2 =- 


6 7 TGL(? 




(13) 


The point on the axis of symmetry of the electron through 
which a perpendicular plane would divide its surface into parts 
having equal charges, will be called the center of charge. If <f) 
is the acceleration of this point, 


Hence 



( 14 ) 


The rest mass m is defined by 


m ~ 


6 7 ra<? 
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Hence, as the dynamical assumption requires that 
K x + K a = 0, 

it follows that the acceleration of the center of charge is deter- 
mined by eE 1 = mf (15) 

at the instant that the electron is at rest in S. 

Consider an electron which has the type of motion under 
discussion, and which is at rest in S' at the instant considered. 
Let the acceleration f of the center of charge make an angle 
with the direction of S” s velocity relative to S. Then, dropping 
the subscript, eE' = mi'. 


Substituting in each of the component equations the values 
of E' x , E' y , E' z from (1), (2), (3), section IQ, and those of/J, 
/»> fl from (27), (28), (29), section 6, it is found that 


eE x = mk 3 f x , 


(16) 

e {A; + i(vxH)„] 

► = mJcf v , 

(17) 

e {^ + i(vxH) 2 ] 

■ = 

(18) 


As the electromagnetic force in S is defined by 
K==e |E + i(vxH)j, 


it is seen to be necessary to distinguish between the longitudinal 


mass 

and the transverse mass 


m l = mJc\ 
m t = mk. 


(19) 

( 20 ) 


Both masses increase with the velocity, becoming infinite as 
the velocity of light is approached. In terms of the transverse 
mass the dynamical equation may be written in the compact form 

K = |o»,o, 

where K is the impressed force. 


( 21 ) 
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21. General case. Consider an electron a point P of which 
is at rest in S at the time 0. Denote by f, f, etc. the acceleration, 
rate of change of acceleration, etc. of this point. Choose axes 
so that the X axis has the direction of f. Then if Q is a neigh- 
boring point of the electron whose coordinates relative to P are 
dx , dy , dz at the time 0, the values of these coordinates at the 
time dt will be 

dx t = dx -f dx + dy + ~ dz) dt 
* \dx dy * dz ) 

and similar expressions for dy t and dz t . But, as / x =/, 

/„= /.= 0 , , 

dz,= dxyjl-^ 

dy, = dy, 

dz t = dz. 


Equating coefficients of like powers of dt and dx, dy, dz in 
the equivalent expressions for dx t , dy t , dz t , it is found that 


^ = ^ = ^ = ^ = = 0 . 

dx dy dz dx 

dx c 2 ’ 


( 22 ) 


*LJLJL 

dy dz dx 


== etc. = 0. 


Hence the velocity is not a function of the coordinates, and when 
one point of the electron is at rest in S, all other points are also 
at rest. Moreover, the y and z components of the acceleration 
are not functions of the coordinates, and the x component is a 
function of x only. 

Before proceeding further, it is convenient to distinguish 
between the orders of possible factors which may appear in the 
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dynamical equation of the electron. If a denotes the radius of 
the electron, it will be considered that 

/3 is of the first order, 

j? 

is of the second order, 
c 

fa 2 

is of the third order, 
c 8 

f a * 

z—r is of the fourth order, and 
c 

is of the fifth order. 

c ? 6 

The dynamical equation to be obtained will be carried through 
the fifth order as thus defined. 

As before, the impressed force on the electron is 

K x = f?E r (23) 

The next step is to evaluate the reaction on the electron of 
its own field. Let the origin be located at a point 0 on the 
surface of the electron, and for the purposes of the following 
analysis let the orientation of the axes relative to f be arbitrary. 
Then if P is another point on the surface of the electron whose 
coordinates relative to O are x, y, z , equation (15), section 16, 
gives for the force exerted by a charge de at 0 on a charge de x at P 

7T ^ dede, f 1 ,3 2 1 3 2 

dK n = 1 < m (1 — v-m H — H — e*m v 

2 * 4 77T 2 V 3 \ 2 y 8 y 8 8 * 

— + - 8*v • . 

15 b ^3 7 / 

- i 7,(i - 1 1 8-m • • •) + 1 8, (i - 2 Y-ra . . •) 

+ (24) 

Integration of this expression with respect to de x will give the 
force exerted on the rest of the electron by the charge de at 0. 
Finally, on integrating with respect to de , the X component of 
the force on the electron due to the reaction of its own field 



52 AN INTRODUCTION TO ELECTRODYNAMICS 


will be obtained. In performing these integrations, the charge 
on the electron may be considered to be uniformly distributed 
over its surface, for, even under the conditions assumed 'in the 
last section, reference to (12) shows that the divergence from 
uniformity there implied leads to no term of less than sixth 
order which does not vanish upon integration. Moreover, it is 
unnecessary to take into account the variation of f from point 
to point of the electron, since (22) shows that the only term 
involved of less than sixth order vanishes when the integration 
is performed. A fortiori the variations of the derivatives of f are 
negligible. 

Omitting from the integrand all terms which vanish on inte- 
gration, (21) leads to the expression 

- i hff** B £(B) + ¥?-¥?(?-$) 

, i f> 7, o**\ 1. 




k 2 =- ? — -J r 

b 7 TdCT 


(* * i 


and similar expressions hold for the y and z components. Hence 
the dynamical assumption leads to the vector equation 
i? e 2 t e 2 • e 2 a •• e 2 a 2 ••• 
eEl = 6W f_ G^ f + 9 (26) 

Consider an electron which is momentarily at rest in S'. 
Then, dropping the subscript, 


-f ' -f --- f ' . 
* 9 77v/ 


f’ .... 


Equations (27), (28), (29), section 6, give 

f x = W x , 

f,j — ky„, 
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Differentiating (24), (25), (26), section 6, once, twice, and 
thrice, with respect to t ', and then placing V r equal to zero, it 
is found that 

f„ = W v + 8W.p^, 


/' = F/ 2 +3/c 5 f-p^; 

f ! x = k b f x + terms of sixth and higher orders, 

f'z = Wf z -\ ; 

/i — • * •> 

“**/;+ •••• 

Substituting in each component of the dynamical equation 
these values of f x , f' y , /', and their derivatives, as well as the 
values of E' x , E' r E[, from (1), (2), (3), section 10, the dynam- 
ical equation of an electron for a system relative to which it 
has a velocity v is found to have the following components : 

,, eVc 3 j e*k 6 

eh x = a o i f ’fi v* • * ’ 

o irac 2 ire 

__ eW u rak* _ eW • v r27 . 

6 8 J J'~' 4 *' * -i o ft J J' * \ J 

7 TO U 7TC lo 7T 6 




e 2 F % m e 2 # 2 & 4 *v 

' 6W> Jv+ Wire* Jv ~ 18 ire 3 1,1 




e 2 /r 2 % ev/F v eVF 
W + 9^T 4 18 7TC 5 


( 29 ) 
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These equations show that it is necessary to distinguish 
between the coefficients of the longitudinal and transverse 
components of the acceleration and of each of its derivatives. 
All the coefficients approach infinity as the velocity of the 
electron approaches that of light, and the series cease to 
converge. 

Including all terms not higher than the third order, the 
dynamical equation of an electron may be put in the form 

K = mf-nf (30) 


where K = e -f ^ (v X H) j- 

is the impressed force, and 


mi 


6 7rac 2 9 


71 = 



22. Rigid body. Consider an element of volume dr of a 
material body large enough to contain a vast number of positive 
and negative electrons, but small compared to the total volume 
of the body. A rigid body will be defined as one all such 
elements of which maintain the same relative configuration and, 
on the average, the same internal constitution with respect to 
the system in which the body is momentarily at rest, whatever 
external conditions the body may be subject to. 

Consider a rigid body momentarily at rest in system S. 
Equations (27), (28), (29) give for the dynamical equation of 
an electron in this body 


e 


{b + 1(txH)| 


f £_ f 

6 rrac 2 6 7T<? 8 ’ 


(31) 


through the third order, which is as far as the analysis will be 
carried in this section. 

The electric and magnetic intensities appearing in this expres- 
sion may be separated into the intensities E e and H e due to the 
impressed field, and the intensities E 0 and H 0 due to the fields 
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of the other electrons in the rigid body. For the latter, equa- 
tions (13) and (14), section 16, give 






where the summation extends over all the electrons except the 
one under consideration. 

Suppose now that there is no external field and that the 
rigid body is permanently at rest in S. Symmetry requires that 
as many electrons in any element of volume dr shall have a 
given velocity or acceleration in one direction due to the internal 
motion as in any other direction. This same condition must be 
satisfied in the presence of an impressed field, for the internal 
constitution of a rigid body is independent, by definition, of 
the external conditions to which it may be subject. 

To return to the case of a rigid body momentarily at rest in 
8 in the presence of an impressed field, let f e and i e be the 
acceleration and rate of change of acceleration of the body as 
a whole, and v 0 , f 0 , and f 0 the velocity, acceleration, and rate of 
change of acceleration of an electron of the body due to the 
internal motion. Summing up over all the electrons, 


X eE o 


i 

6 7T6* 2 



iVev o XH = 0, 

c ** 


1 

6 7 rc* 



and therefore 


<**> 

where the double summation is for all values of both i and j 
such that / 7 4 /3 
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Consider a rigid body at rest in S'. Then, dropping the 
subscript e, 

_ { S' 1 ’ + ' " 


Substituting for E f its value in terms of E and H, and for f f 
and f f their values in terms of f and f, the dynamical equation 
of a rigid body takes the form 



The first term in the brace multiplying f is the sum of the 
masses of the electrons composing the rigid body, while the 
second term is the sum of the mutual masses of these electrons. 
While the mass of an electron must be positive, the mutual 
mass of two may be positive or negative according as they have 
like or unlike signs. Hence the mass of a rigid body is greater 
the more electrons there are of the same sign. The same is true 
of the coefficient of the rate of change of acceleration. In fact 
this coefficient vanishes if the body is uncharged. 

Denoting by K the impressed force, and by m and n the con- 
stant coefficients of f and f respectively, the dynamical equation 
of a rigid body takes the form 

K = 77lf — 7lf • • •• (34) 


A conductor carrying a current may be considered, in so far 
as the expressions arrived at in this section are concerned, as one 
rigid body through which another is passing. As the electrons 
carrying the current are all of the same sign, their mutual masses 
are positive, and the mass of the current is greater than the sum 
of the masses of the individual electrons which constitute it. 

23. Experimental determination of charge and mass of electron. 
Consider a stream of negative electrons from the cathode of a 
discharge tube travelling at right angles to the lines of force of 
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mutually perpendicular electric and magnetic fields, 
magnetic force on each electron is given by 


Hence if 


K = 


-v X H = 


e |e + -v 
-E, 


c 


X 



The electro- 


this force will vanish, and 


E 

v = — c. 
1L 


By adjusting crossed electric and magnetic fields so as to 
produce no deflection in a beam of cathode rays, J. J. Thomson 
has found the velocity of these charged particles to be about 
one-tenth that of light. 

If, now, the magnetic field be suppressed, 

K = eE, 

and to a first approximation 

mi = eE. 


If the rays suffer a deflection d in travelling a distance 8 
through this field, 


, 2 dv* 

„2 ’ 


and the ratio of charge to mass is given in terms of measurable 
quantities by the expression 

e 2 dv 3 


m Es 2 


(35) 


Similarly, if the electric field be suppressed, the value of the 
ratio of charge to mass may be obtained from the deflection 
suffered in traversing a magnetic field. In this case 


e __ 2 dvc 
m Ms 2 


(36) 


By these methods it is found that 

- = 1.77(10y<?V4w (37) 

m 

for the negative electron. 
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Determinations of this ratio for beta rays moving with veloci- 
ties only slightly less than the velocity of light verify the theo- 
retical expression (20) for the increase of transverse mass with 
velocity. 

In order to measure e and hence m other experimental 
methods are necessary. Suppose an electron to be attached to 
a minute oil drop situated between the horizontal plates of a par- 
allel plate condenser. If the electric field is adjusted so that the 
oil drop remains at rest, its weight w is balanced by the force 
eE. Now if the drop is allowed to fall freely through the sur- 
rounding gas, its radius may be calculated from its rate of fall 
by Stokes’ law. From the radius and density of the oil w may 
be determined and hence e computed. In this way the electronic 
charge has been found by Millikan to be 

e = 4.77 (10) -10 V47r, (38) 


and, combining this with (37), the mass of the negative electron 


is found to be 


m — 9.0 (10) -28 gm., 


(39) 


which is about one eighteen-hundredth of the mass of a hydro- 
gen atom. Hence the radius of the negative electron is 


a = 1.88 (10)~ 18 cm. 


Since the positive electron has not yet been isolated, it has 
been impossible to measure its mass and radius, although there 
are reasons for supposing that it has a much greater mass and 
consequently a much smaller radius than the negative electron. 



CHAPTER Y 


EQUATIONS OF THE ELECTROMAGNETIC FIELD 

24. Divergence equations. The field due to a point charge is 
completely specified by any two of the three vector functions of 
position in space and time, E, H, and c. In the case of a com- 
plex field, E and H are the resultants of the electric and mag- 
netic intensities respectively of the component simple fields, but 
c must be given for each elementary field. Hence in order to 
avoid explicit reference to the components of a complex field, 
as well as in order to give the field equations as great a sym- 
metry as possible, the field is usually described in terms of 
E and H. 

To find the divergence of the electric intensity due to any 
given distribution of point charges consider a small region dr 
surrounded by the closed surface a. By Gauss’ theorem, 

V*E dr— f E •da. 


Let letters with strokes over them refer to component fields. 


Then E d<T =~ 

for each element of charge. 
Therefore V • E dr = 


dN 

da 2 


d<T 



Now the part of this sum due to charges outside the region 
dr vanishes, while the part due to the charge de inside this 
region becomes equal to de itself. Hence 

V-E = p, (1) 

where p is the density of charge at the point in question. 

• 69 
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In section 13 it was shown that 

H=V x a. 

Therefore V*H = 0 (2) 

identically, or the divergence of the magnetic intensity is every- 
where zero. 

25. Vector fields. Any vector function of position in space 
and time may be represented by moving lines such as have 
been employed to give a geometrical significance to the electric 
intensity. These lines will be continuous at all points where 
the divergence of the vector function vanishes, as is obvious 
from the discussion contained in the preceding section. At 
other points lines will either begin or end. Equation (1) shows 
that the resultant electric intensity may be represented by con- 
tinuous lines at all points except those at which electricity is 
present, while equation (2) shows that the resultant magnetic 
intensity may be represented everywhere by continuous lines. 

Let V be a vector function whose magnitude and direction 
are represented by lines all points of which are moving with 
velocities of the same magnitude c. If dN tubes of these lines 
(a tube being a bundle of M lines) pass through a small surface 
ds x with normal parallel to the X axis, 

* A, 


For the moment assume that no new lines are formed. Then 
in a time dt , V x may suffer a change due to three causes. In 
the first place, the number of lines passing through ds x may 
increase by virtue of the fact that the lines whose motion will 
bring them through this surface at the end of the time dt are 
more closely packed than those passing through ds x originally. 
The increase in V x due to this cause is 



dV x , dV u 

— -c H 3 

dx x dy 


8 J* C + d l* e \ 

dv y dz Z S 


dt . 


Secondly, the velocity associated with the new lines may have 
a different direction from that of the old. This will produce a 
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crowding of the lines during the time dt and account for a 
change in V x equal to 


-{ 


rr d C » . -rr 

r -^ +F -&i 


dt. 


Finally, if c differs in direction at neighboring points on the 
same line, there will ensue a twisting of the lines which will 
produce in V x an increase 



dt 


in the time dt. 

Therefore the total rate of change of V x is 


dt 


c-vr.-F.V-c+V-V*,, 


or 


dV 

— -f cV*V == V x {c x V}. 

ct 


( 8 ) 


Now consider the increase in V due to the formation of 
new lines. Attention will be confined to those fields whose 
lines terminate only on sources. Let 
the points O etc. in Fig. 9 be each the 
source of a new line emitted in the direc- 
tion of the arrows, the line sources them- 
selves having a velocity v to the right. 

The number of tube sources per unit 
volume is obviously 

V-V, 

and as the direction of the lines at P 
is QP, the increase in the value of 

V at P in a time dt due to the formation of new lines is 

(c — v)V*Veft. 

Therefore the complete expression for the rate of change of 

V becomes 

|L + vV-V=Vx{cxV}. (4) 

ct 
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If the field under consideration is due to the superposition 
of a number of simple fields of the character of those just dis- 
cussed, ™ 

+2vV-V= V x %{c x V>, (5) 

where the stroked letters refer to the component fields. More- 
over, since a number of point sources emitting lines in different 
directions may be considered to coalesce, this equation applies 
equally well to point sources from which lines diverge in all 
directions. It is to be noted that equation (5) is a consequence 
of the properties of three-dimensional space, and nothing more. 

26. Curl equations. Replace V by E in (5). Substituting 
from (1), the second term of the left-hand member becomes 

'Z'pv- 

Since, however, only one element of charge may occupy one 
point in space at a given time, the summation sign may be 
dropped. 

The right-hand side of (5) becomes 
Vx^{c £ xE>, 

where c E is the velocity of a moving element of an electric line 
of force of one of the component fields. 

But H = - ^ c £ x E. 

Hence VxH = -{E + pv}. (6) 

c 


If V is replaced by H in (5), it follows from (2) that the 
second term of the left-hand member vanishes. Hence the 
equation becomes 

* -L r — ^ x w (j?') 


VxSi{5*xH}=lH, 


where the velocity of a moving element of a magnetic line 
of force of one of the component fields does not in general havQ 
the same direction as the velocity c E of a moving element of an 
electric line of force of the same field. 



EQUATIONS OF ELECTROMAGNETIC FIELD 63 


But if a is eliminated from the equations obtained by curling 
(20) and differentiating (21) with respect to the time, in section 18, 

it is found that i # 

V X E = — - H. (8) 


Hence, as (7) and (8) apply equally well to simple or com- 
plex fields, comparison of these two equations shows that for a 

single elementary field _ i _ _ 

E=--5 ff xH+Vf, (9) 

a relation which is complementary to the definition 

E. (10) 

Therefore E=-^c fl x(c £ xE)+ 

or l(i- =- *»+▼+ s 

and H = — & E X (Cj^- x H) + — c E x V 

or h(l-££^ = 1c £ x Vf. 

Hence it follows that % is parallel to c E when, and only 
when, V^r is parallel to c E . 

27. Electrodynamic equations. Equations (1), (2), (6), (8) 
specify the electromagnetic field in terms of the distribution 
of density of charge p and velocity v. If (2), section 18, is 
written as an equation instead of as a definition in order to 
signify that it includes the dynamical assumption immediately 
following it, this equation suffices to determine the effect of 
an electromagnetic field upon matter. These five equations 
contain the whole of electrodynamics. Collected they are: 

V-E = />, (11) V-H = 0, (18) 

V x E = --H, (12) VxH = l{E + pv}, (14) 

G Q 

F = p{e + 1(vxH)|» ‘ (15) 

where F is the electromagnetic force per unit volume. 
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28. Energy relations. From (12) and (14) it follows that 

H*V X E — E*V X H = — - (E-E + H-H)— - pE-v. 

0 c 

But H*V x E — E*V x H = V*(E x H), 

and pE*v=F»v. 

Hence ^ {- (A 2 + -ff 4 )} + «V-(E x H) + F.v = 0. 


Integrating over any arbitrarily chosen portion of space, and 
applying Gauss’ theorem to the second term, 


Jtf\( Ei + H ^ dT + C f ( ExH W <r+ f F*vdr = 0, (16) 


where dr is an element of volume, and d<r an element of the 
bounding surface having the direction of the outward drawn 
normal. The third term of this expression measures the rate at 
which work is done by the electromagnetic field on the matter 
contained in the region selected. Hence, if the principle of con- 
servation of energy is to hold, the first two terms must be 
interpreted as the rate at which the energy of the field increases 
plus the rate at which energy escapes through the surface 
bounding the field. 

Suppose E and H to be zero everywhere on this surface. 
Then no energy escapes from the region enclosed, and the rate 
at which work is done on the electromagnetic field equals the 
rate at which its energy increases. But the integrand of the 
second term of (16) is everywhere zero. Hence the rate of 
increase of energy of the field must be represented by the first 
term, the form of which suggests that 


u = $(E* + H*) 

is to be considered as the electromagnetic energy per unit volume. 

The second term of (16), then, must be interpreted as the out- 
ward flux of energy through the surface enveloping the field, and 
its form suggests that rF m 


is to be considered as the flow of energy per unit cross section 
per unit time. 
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29. Electromagnetic waves in space. For empty space the 
electromagnetic equations take the form 

V-E = 0, (17) V*H = 0, (19) 

VxE = -1h, (18) V X H = i E. (20) 

c c 

To eliminate H, curl (18) and differentiate (20) with respect 
to the time. Thus 

1 

V x V x E = V x H 

c 



But Vx V x E ~ VV-E — V«VE 



= - V-VE. 


Therefore 

V*VE — ijj = 0, 

c 2 

(21) 

and, similarly, 

V-VH — = 

c 2 

(22) 


These are equations of waves moving with velocity c. 
Consider a plane wave advancing along the X axis. Then E 
is a function of x and t only, and (21) reduces to 


a 2 E l a 2 E 

dx 2 (? dt 2 


(23) 


The solution of this equation for the case of a wave moving 
in the positive X direction is 

E x 

E z = h(x — cf). 

Hence it follows from (17) that E x is a constant, and since 
the present discussion is concerned only with the variable part 
of the field, this constant may be taken as zero. 
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Therefore Ji x = 0, (24) 

E,= g(x-cf), (25) 

E z = h (x — ct'). (26) 

From (18) it follows that 

0 , 

H y — eh r (x — cf), 

H = -cg'(x-c€), 

and hence H x — 0, (27) 

Jf y = — h(x — c tf), (28) 

H=g(x-ct'), (29) 

except for a possible constant of integration. Therefore the 
variable parts of E and II have the same magnitude at any 
point and time, and lie in a plane at right angles to the direc- 
tion of propagation. 

The cosine of the angle between E and H is proportional to 
E y II y +E z II = -gh+gh 
= 0 , 

showing that E and H are at right angles in the wave front. 
The preceding section gave the flow of energy as 

c(ExH)=</+y)i 

= c {^(^ 2+7/2 )} i ’ (30) 

showing that the propagation of energy is along the X axis and 
that the entire energy of the wave front is advancing with 
the velocity of light. It follows from this equation that the 
direction of propagation of the wave is at right angles to the 
plane of E and H in the sense of 


E x H. 
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Consider a spherical wave having its center at the origin. In 
this case E is a function of the radius vector r and t only, 
and (21) reduces to 


2aE_l^E 

dr 2 r dr c 2 dt 2 


(31) 


of which the solution is 

E~- r f(r±ct), 


E e = ~ff(.r± ef), 
r 

E <> —^-h(r± ct~). 


As in the case of the plane wave, it may be shown that 
E and H are mutually perpendicular, and at right angles to the 
direction of propagation. If r 1 is a unit vector along the out- 
ward drawn normal to the wave front, the flow of energy is 

(T 2 _L Jl 2 

e (E x H)= ± c g -±± r. 


= ±c-{|(7? 2 +i/ 2 )}r 1 , (32) 

showing that the amount of energy passing through unit cross 
section in unit time varies inversely with the square of the dis- 
tance from the source, and that the entire energy of the wave 
front is advancing with the velocity of light. 

30. Radiation pressure. Substituting in (15) the values of p 
and pv from (11) and (14), the electromagnetic force per unit 
volume takes the form 

1 • 

F = V*EE + (V x H) x H E x H, 

c 

and making use of (12) and (13), 

F = — 1 — (ExH)+V*EE+V*HH+(VxE)xE+(V xH)xH. 

c dt 
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The total electromagnetic force on the matter in any given 

region is r 

K = I Fdr 

= 1^+ K ff , 


where 


1 dr 
<? dtj S 


and 




sdr, (33) 

U U IJ 

= J {V*EE+V-HH + (VxE)xE+(VxH)xH}dr. 
Now let J = J { V-EE + (V X E) x E} dr 

= j {V'EE + E-VE - i V(E*E)}c?t 

=j (V-(EE)- V(E-E)} dr. 

By Gauss’ theorem, 

J V*(EE)cZt = J EE -d<r, 


and 


J V (E*E) dr =i j V*(iE-E)<ir + ete. 
— i J E*Ei*i 


•c?cr + ete. 
Edar. 


=/e- 

Therefore J =J* EEvfo* — | J^E^Edcr, 
and K a is given entirely by the surface integral 

(EE + HH)*cZ<r — | j (E«E + H*H)cZ<r (34) 

taken over the surface bounding the chosen region. 

Consider a region which is the seat of either stationary radia- 
tion or undamped periodic vibrations. In the former case the 
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average K T is zero, and in the latter case K T vanishes provided 
the value of the force desired is the average over a whole num- 
ber of periods. Hence the force on the matter within the region 
is given entirely by K a . 

Consider a box (Fig. 10) with perfectly reflecting walls con- 
taining homogeneous isotropic radiation. Describe the pill-box 
shaped surface ABCD about an element da of that part of the 
wall of the box which is perpendicular to the X axis. The only 



matter within this closed surface is the portion MN of the wall 
of the box. In so far as the effect of the radiation inside the 
box is concerned, the average force on MN is given entirely by 
K a , which reduces to an integral over BC, since AD lies com- 
pletely outside the field, and the sides AB, CD of the pill-box 
have a negligible area. Hence if X x , Y x , Z x are the forces per 
unit area on the surface MN of the box due to the pressure of 
the radiation inside it, parallel respectively to the A, F, Z axes, 

Y=E x E y +HJI y , (36) 

Z=E x E z +H x H z , (37) 

and similar expressions hold for the parts of the wall perpen- 
dicular to the Y and Z axes. 
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As the radiation is assumed homogeneous and isotropic, 

E X = E=E S = 0 , 

where the stroke denotes the average value of the quantity over 
which it is placed. Similar relations hold for the components 
of the magnetic intensity. 

Therefore 

x x =- 1 (¥>+!?-) 

Y X =Z X = 0, (39) 

showing that the radiation 
within the box exerts a 
normal stress on the walls 
equal to one third of the 
electromagnetic energy per 
unit volume. 

Consider, now, a box 
similar in shape to that of 
Fig. 10 with perfectly con- 
ducting walls. Instead of being filled with homogeneous radia- 
tion, suppose the box to contain a train of undamped plane 
waves which are incident on the surface MR (Fig. 11) at the 
angle 6 . 

Referred to the X l Y 1 Z 1 axes in the figure, equations (24) to 
(29) give for the incident wave 

K = o, 

\= - ct )> 

\(X - ct), 

K = o, 

II z=r/X x - <;t ')- 



Fig. 11 


and 
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Therefore the components of the electric and magnetic 
intensities along the XYZ axes due to the incident wave 
alone, are E= 9l mx6, 

Ey=-g i WSd, 

If x = — h l sin 0, 

H y — h x cos0, 

Hz— 9 v 

where the argument of the arbitrary g l and Ti x functions is 


x. 


ct = — x cos 6 — y sin 6 — at. 


Similarly, for the reflected wave, 

S x -g t sai6 % 

Ey = g 2 GOS0, 

U x = — \ sin 0, 

H v = - \ cos 0 , 

H z = //o, 

where the argument of the g 2 and h 2 functions is 
x 2 — ct = x cos 0 — y sin 0 — ct. 

At the reflecting surface the arguments of the functions with 
different subscripts become the same. Moreover, as this surface 
is perfectly conducting, the tangential component of the resultant 
electric intensity must vanish. Therefore 

— g x cos 0 9 2 cos 0 = 0, 

^i+A 2 =0, 


or 


9 1 ~ 9 2 = 9 > 


Hence it follows from (30) that all the energy brought up 
in the incident wave is carried away in the reflected wave. 
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Therefore a perfectly conducting surface is a perfect reflector of 
electromagnetic waves. 

The resultant electric and magnetic intensities just to the 


right of MN are 74= 2 <7 sin 0, (40) 

^=0, (41) 

a;= 0, (42) 

ir,= 0, (43) 

H y = 2 7t cos 0, (44) 

7A= 2//. (45) 


Hence, as there is no field to the left of MN, this surface 
must have a charge 2 g sin 6 per unit area. If g is a simple 
harmonic function, this charge will be alternately positive and 
negative. 

Substituting the values (40) to (45) in the expressions (35), 
(36), (37) for the components of the stress on the reflecting 

surface ’ X x =- 2 (/+ A 2 ) cos 2 0, (46) 

= (47) 

showing that the stress due to the radiation is a pressure nor- 
mal to the surface. 

Now, the energy density of the incident wave is 

and therefore X x =— 2 u x cos 2 6. (48) 

Let / stand for the energy striking the reflecting surface per 
unit area per unit time. Then 

I =u t c cos 0, 

and X x —~ 2 ^cos0. (49) 

Suppose that instead of containing a single train of plane 
waves, the box is filled with plane waves travelling in all direc- 
tions at random ; that is, with a homogeneous isotropic radiation. 
Then if u is the energy density of the radiation, the energy per 
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unit volume of that portion of the radiation which is incident 
on MN between the angles 6 and 0+ d0 is 

w x = sin#d0, 

and X*= — 2^)w 1 cos 2 0 


= -‘jf 


cos 2 0 sin Odd 


1 

3 U ' 

agreeing with (38). 

31 . Electromagnetic momentum. Consider a closed surface 
ABCB (Fig. 12), surrounding some matter. Let the field inside 
this surface be in a stationary state. Then the force on the matter 
contained is given by integrated over the bounding surface, 
the outward-drawn 
normal being posi- 
tive. Denote by 
K al the part of K a 
obtained by inte- 
grating over B C. 

Then the force K 
on the matter under 
consideration, due 
to the electromag- 
netic field extending 
through the surface 
BC, is equal to K^. 

Let BCEF be a closed surface which surrounds no matter, 
such that the radiation field inside it comes into contact with the 
surface only between B and C. Then the value of integrated 
over this surface is equal to since the outward-drawn 

normal to BC in this case has the opposite direction to that for 
ABCB. Hence __ jr k — o, 

as no matter is present inside the surface. Therefore 

k-k t =o, 



Fig. 12 


(50) 
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where K is the force exerted by the electromagnetic field in 
BGEF on the matter in ABGD. Hence, if the law of action and 
reaction is to hold, — K T must be interpreted as the force exerted 
by the matter in ABCD on the electromagnetic field in BGEF. 
The form of this expression, namely, 


suggests that 



1 


s 


& 


is to be considered as the momentum per unit volume of the 
electromagnetic field. 

Suppose the matter in ABC!) to consist of a rigid body with 
a plane, perfectly reflecting surface just to the left of BC. Let 
a limited train of plane waves be incident on this surface at 
the angle 6. Reference to (30) shows that the electromagnetic 
momentum per unit volume of the incident radiation is 



c 


( 51 ) 


Therefore the momentum of the radiation striking each unit 
area of the surface in unit time is u x cos (9, and the momentum 
of the reflected radiation is of the same magnitude. The force 
exerted by each unit area of the reflecting surface on the train 
of waves is equal to the vector increase in momentum per unit 
time ; that is, — K r = 2 u x cos 2 0 

along the outward-drawn normal. Consequently the stress 
exerted by the radiation on this surface is 

K = — 2?qcos 2 0 

= — 2 -cos#, 
c 

agreeing with (48) and (49). 

32. Four-dimensional representation. In four-dimensional space 
two mutually perpendicular lines may be drawn at right angles 
to an element of surface. Consequently the vector properties 
of such an element cannot be expressed by the direction of its 
normal. In fact, it is necessary to distinguish between directed 



EQUATIONS OF ELECTROMAGNETIC FIELD 75 


linear elements, directed surface elements, and directed volume 
elements. The first have as components their projections on 
the four coordinate axes, the second their projections on the 
six coordinate planes, and the third their projections on the four 
coordinate planoids. Hence the first and third are often called 
four-vectors, and the second six-vectors. Here they will be called 
vectors of the first, second, and third orders respectively. 

Let a, b, and c be three vectors of the first order. Then the 
vector product a X b is a directed surface having for its area the 
parallelogram of which a and b are the sides, and so directed that 

bxa = — axb. 


Similarly, a x b x c is a directed volume having for its magni- 
tude that of the parallelepiped whose edges are a, b, and c. In 
the case of any cross product, interchange of two adjacent vectors 
changes the sign of the product. 

Lot k x , k 2 , k 3 , k 4 be unit vectors of the first order parallel 
respectively to the AT, F, F, L axes. Then 

^23 “ k 2 x ^3 

is a unit vector of the second order in the YZ coordinate plane, and 

k m = k S Xk 4 Xk l 


= k 34 Xk ! 


is a unit vector of the third order in the ZLX coordinate 
planoid. From their definition it is obvious that two adjacent 
subscripts of a unit vector of the second or third order may be 
interchanged provided the sign of the vector is changed. 

The dot product of two unit vectors is defined as follows. 
If the vector of lower order has a subscript which the other 
lacks, the dot product vanishes. Otherwise, the product is the 
unit vector remaining after like digits in the subscript of each 
vector have been brought to the end and cancelled. Thus 


^ 124*^12 ^ 412*^12 
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Thus, while the order of the cross product of two unit vectors 
is equal to the sum of their orders, the order of the dot product 
is equal to the difference of the orders of the two factors. 
Consider the vector of the second order 

m = m jv \, + m „*K 3 + + M *iK + M «iK * + ( 52 ) 


The dual M x of M is defined as a vector of the same order 
with components such that 

where mnop is formed from xyzl by an even number of inter- 
changes of adjacent letters. Hence 

M* = M a . k 12 + MJ k J3 + M y , k 31 + M az k l4 + 3/„k 24 + M n , k, 4 . (58) 
If p = iik t + r„ k 2 + Jik 3 + 1\ k 4 , (54) 

the rule for forming the dot product shows that 

P-M = ( +llM x -P z M m ,+ J’MK 

+ (- P X M XII + i;m vz + P t M yl ) k 2 

+ ~ P „ M yz +«)k 3 

+ (- P x M xl - P y M yl - P z M zl )k 4 . (55) 


In four-dimensional analysis the vector operator 



plays much the same part as V in three dimensions. 
Consider the product 


= / 








r „ 

oy 

dz 

+ dl ) 

, / 


1 

dM m 

yz 


\k 

+ \ 

K dx 

1 

dz 

dl j 

r i 


' ZKr 

K dx 

dM az 

dy 


dl t 

) k 3 

. I 





w 

+ \ 

\ dx 

~dy 

dz 


) 4 


( 56 ) 
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Remembering that l = ict, equations (11) and (14) of the 
electromagnetic field may be written 


oil 8H„ 8iE x 1 

+ ~xf = ~pv x , 

8y 8z 81 c 




d JL_ 


1 

8x 

811 v 

i 

8H X 

dz 

81 

8iE z 

~e P * 

1 

dx 

8iE x 

dy 

8iE 

diE z 

81 

= -pv z , 

= ip, 

dx 

' Q 

dy 

dz 



and equations (12) and (13) 


8iE\ 

dx 



d JJL± d JEz+. d JL 

8y^8z 81 

djK 8H U 

8z 81 

8iE x 8H r 

8y 81 


= 0, 
= 0, 
= 0, 


811 8H„ 8H, 

? IL — ? 

dx dy dz 


Comparison with (52), (53), (54), and (56) shows that if 
M ee H z k 12 + H x k 23 + 11 y k 31 - iEX* “ iE vK* ~ 

and P e= i ^ p^k 2 + - pv z k 3 + «>k 4 , 

the two scalar equations (11) and (13) and the two vector equa- 
tions (12) and (14) of the electromagnetic field are expressed 
by the pair of four-dimensional vector equations 


0*M = P, 

(57) 

o 

• 

K 

X 

II 

p 

(58) 
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The equation (15), which gives the effect of an electro- 
magnetic field on matter, and the energy equation may be written 
together in the form 

F x = +\pv y Ih-\pvJI v+9 E x , 

F„=- - pv x H z +~pvj' x + pF y , 

6 6 

*\= ]pv x H y -\pv y H x + pE„ 

i du 1 , 1 . „ , 1 . T , 

-cdt = -P^F x +-pvsE y +-pv t ih,. 

Hence, if F = F^-t- F y \+ F t \+ l - ^k 4 , 

these relations are contained in the four-dimensional vector 
eq uafci °n F = P-M. (59) 

The scalar operator 

A <?_ d 2 _ 

+ dy l + az 2 " l " dl 2 

is known as the d’Alembertian. In four-dimensional analysis 
the wave equations (21) and (22) are expressed by 

0-0M = 0. (60) 

In addition to those just mentioned, many other electrody- 
liainic relations may be expressed in far more compact form in 
terms of four-dimensional vector analysis than in the analysis of 
three dimensions. 



CHAPTER VI 


RADIATION 


33. Radiation from a single electron. Let v and f be the 

electron’s velocity and acceleration respectively at the time t. 
Take its position at this time as origin, and choose axes so that 
the X axis is parallel to v and the XY plane contains f. Describe 
about the origin a sphere of radius r, large compared with the 
radius of the electron. Consider an element d<r of the surface 
of this sphere such that the radius vector drawn to it from the 
origin makes an angle 6 with the X axis. The energy emitted 
from the electron during an interval of time dt in the direction 
of this radius vector will reach the surface of the sphere at a 
time t + r/c, and will take a time 

(1 — /3 cos 0) dt 

to pass through this surface. Hence, as the flow of energy per 
unit cross section per unit time is equal to 
s=(?(ExH), 

the energy emitted by the electron during the time dt is given 

by the integral n 

Rdt == / s*dc (1 — /3 cos 6 ) dt 

taken over the surface of the sphere. Therefore the rate of 
radiation from the electron is 


As 

it follows that 


R =J* &»dcr (1 — /Scos 0 s ). 


H == - c x E,* 

c 

s = E* c — E»cE. 


( 1 ) 


The sphere over which the surface integral is to be taken 
may be made as large as desired. Consequently the terms in 
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equation (14), section 12, for the electric intensity which involve 
the inverse square of the radius vector may be made so small 
compared with the term involving the inverse first power of 
this, quantity that they may be neglected. As the term involv- 
ing the inverse first power defines a vector perpendicular to c, 


s = J£ 2 c 


e 2 C f 


r 

, 2 ftfj % 

a-/n/? v 

(1— j8cos0) 4 

(1— /3cos0) 6 

(l-iScos^yi 


where f r is the component of the electron’s acceleration along the 
radius vector. For small values of /3 this becomes approximately 


e 2 / 2 sin 2 a 
16 7rW 


where a is the angle which the radius vector makes with f. 
Hence the radiation vanishes in the direction of f and is maxi- 
mum at right angles to this direction. 

To find the rate of total radiation from the electron, integral 
(1) must be evaluated. Substituting the expression for s given 

by (2), 


B 16 7rV {^/r 2 (l-/3cos0) a + 2/3 f r\ 1- 


fJM 


ft cos#) 4 


e 2 




fa? 


fttll 


— (1 — # 2 ) f \ 

K V r*(l-ftcoseyj 


6 ire* \(l-fty (1-/3 2 ) 


} 


(3) 


If f' is the acceleration of the electron at the time t relative 
to the system in which it is, at that instant, momentarily at 
rest, reference to equations (27), (28), (29), section 6, shows that 


B 


e*f' 2 

6 77 m C* 


(4) 


34. Radiation from a group of electrons. In the calculations 
of this section it will be assumed that the greatest distance 
between two electrons of the group is small compared with the 
wave length of the radiation emitted, and that the velocity of 
the fastest moving electron of the group is small compared with 
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that of light. If a sphere of radius r, large compared with the 
dimensions of the group, is described about some point in the 
group as center, the energy passing through the surface of this 
sphere in unit time is found as in the last article to be 





( 5 ) 


where 


s = E\ 


= cgE,-E j, 

V 

the summation being taken over the n electrons in the group. 

Now = ejfi-cc-e^} Q 

' 4 vrc* K y 


at a great distance from the ith electron. 

Therefore 8 = lei^V S f < ,cf / c ) • ( 7 ) 


Suppose that the sum of the components of the accelerations 
in any direction is equal to that in the opposite direction. 
Then s vanishes for all directions of c. Hence a ring of any 
number of evenly spaced electrons which are rotating about a 
common axis with constant speed, will emit no appreciable 
radiation. 

To get the total radiation from the group of electrons, sub- 
stitute (7) in (5) and integrate. In this way it is found that 


E 



( 8 ) 


35. Energy of a moving electron. The dynamical equation of 
an electron has been shown to be 

K = mi — nf, 


through terms of the third order, where 

e 2 - 
W ~6t ra<?' 
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The work done on the electron in a time t by the impressed 
field is 

W = J K-vdt 

= m j* t'Vdt — n J" bvdt 

= \ m « - V D - n (y-A ~ v r f .) + f f* dt > ( 9 ) 


where the subscript 1 denotes the value of the quantity to 
which it is attached at the beginning of the time t and the sub- 
script 2 the value of this quantity at the end of this time. 
The first term in this expression represents the kinetic energy 
of the electron, the second term its acceleration energy, and the 
third term the energy which has been radiated. The first two 
terms are reversible in the sense that the energy which they 
represent may be recovered when the electron returns to its 
original state of motion, but the third is irreversible. 

Consider an electron which starts from rest and acquires a 
velocity v by virtue of a very small acceleration continued for 
a very long time. By making the derivatives of the accelera- 
tion small enough, the second and succeeding terms of the 
dynamical equation may be made as small as desired com- 
pared with the first term. Hence, if account is taken of the 
variation of mass with velocity, this equation is given to any 
desired degree of accuracy by 


K = 


m 


r* 


( 1-/3 2 )* 

and the work done is expressed as closely as may be desired by 


W=j‘K-vdt 

■ m<? 


i i .a 

J 


6 7m[Vl— /3 2 


-i • 


( 10 ) 
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In this case the energy radiated is inappreciable, so all 
the energy acquired by the electron may be recovered if 
it is brought back to its original state of rest by a similar 
process. 

Now the energy of an electron moving with constant velocity 
v is given by 

U =\f (E 2 + 11*) dr. 


Inside the electron’s surface both E and H vanish, while out- 
side the surface the values of these intensities are given by 
equations (3) and (4), section 19. Hence if the angle which the 
radius vector makes with the direction of the electron’s velocity 
is denoted by 0, 


i 2 (l — /3 2 ) 2 rr l + /3 2 sin 2 fl fl . 

16 7r JJ r\ l-/3 2 si: 

where r goes trom - = ■ = = 

6 Vl— /3 2 sin 2 # 

to infinity, and 0 from 0 to tt. 

Integrating, U = ^{^== - J 

and the increase in energy is given by 

TT — TT — — — I ivi — /3 2 — ?1. 

0 6 7r4vr=^ 4 4J 


• 2 ^s s ' m0 d0dr, 

sin By 


( 11 ) 


( 12 ) 


The discrepancy between equations (10) and (12) arises from 
the fact that in the calculation from which the former was 
obtained no account was taken of the work done against the 
electron’s field in connection with the progressive contraction 
which takes place as its velocity increases. In order to deter- 
mine the work done in this process of contraction, it is neces- 
sary to evaluate the stress K on each unit area of the electron’s 
surface. From the expressions for the electric and magnetic 
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intensities given in section 19 it follows that the electromag- 
netic force per unit volume just outside the surface is given by 


<5. 

II 

: 

r 


pe (1— -/3 2 )eos# 

_ 4 7rr 2 (1-/3W#/ 


A 


F = p{E-$H} 

/ 

yf\ 


pe (1 — J3*y sin 6 

1 / 


i 

“ 4?rr 2 (l-^sin 2 #) 1 ’ 

' i 

A \ 

L, 

where the X axis is taken in the 
direction of the electron’s velocity. 

If a (Fig. 13) is the angle which 
the normal to the surface makes 
with the X axis, 


J 


tan a = (1 — /3 2 ) tan 6. Fig. 13 


Therefore the tangential component F t of F is given by 
F t = — F x sin a -f- F y cos a 


= 0, 

and the normal component F n by 
F n = F r cos a +F V sin a 


__ pe ( 1 — /3 2 ) V cos 2 4- ( 1 —• yQ 2 ) 2 sin 2 # 
“4 7rr 2 (l-/3 2 sin 2 6>) 4 


(13) 


(14) 


Now the charge per unit area is easily shown to be 

« ( 1 -/ 3 2 ) 


4 ?rr 2 Vl — /S 2 sin 2 0 V cos 2 6 + (1 — /8 2 ) 2 sin 2 6 


( 15 ) 


Therefore, as the electromagnetic force just inside the electron’s 
surface is zero, 

K= 2 ^;+°) 


e 2 (1-/S 2 ) 2 
32 ir 2 r* (1 — ;8 2 sin 2 #) 2 

2 

e 

32-TrV' 


(16) 
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As this stress has the nature of a hydrostatic tension inde- 
pendent of the velocity of the electron, Poincare has been led 
to suggest that the electron may be held together by an equal 
and opposite hydrostatic pressure of an extra-electrodynamical 
character. 

As the electron’s velocity increases from 0 to v , its volume 
decreases from 4 » 

i 7 ™ 

to |tiW1-/3 2 . 

Therefore the work done against the stress of (16) is 

< i? > 

Adding this to (10), (12) is obtained. 

36. Diffraction of X Rays. Equation (3) of section 33 shows 
that in general an electron will radiate energy whenever it is 
accelerated. However irregular its motion may be, the radiation 
emitted may be analyzed by Fourier’s method into a series 
of superimposed simple harmonic waves. Waves of a length 
from 4000 A to 8000 A constitute light of the visible spectrum, 
whereas waves of a length of the order of 1 A are called X rays. 
These rays have great penetrating power, and all attempts to 
diffract them were unsuccessful until Laue suggested in 1913 
that the distances between adjacent atoms in crystals were of 
such a magnitude as to make these substances suitable natural 
gratings for the diffraction of X rays. The following theory is 
presented very nearly in the form given originally by Laue. 

Let a r a 2 , a 3 be vectors having the lengths and directions 
of the edges of an elementary parallelepiped of the crystal. 
Then if x, y , z are the coordinates of an atom relative to an 
origin 0 at the center of the crystal, 

X = ma lx +na 2x +pa zx , 

y = ma ly +na iy +pa ail , (18) 

Z = ma lx +na. 2z +pa tz , 

where m , n, p are positive or negative integers. Let r be the 
distance of this atom from the observer at -P, and B the distance 
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of the center of the crystal 0 from P. These two distances 
will, in general, be so nearly equal that they may be considered 
the same, except in so far as the phase of the radiation is 
concerned. 

Let the incident radiation be plane, the direction cosines of 
the wave normal being denoted by a Q , /3 0 , y Q . Then if the intra- 
atomic vibrators are all alike, the displacement of any vibrator 
at a time t will be given by the real part of 

(x« 0 + Vf3o + *yo- ct ) 

A 0 e A 

and the field strength at P due to this displacement by 

A (r + + vfio + Wo “ c 0 
-e * 


since the electric intensity at a distance from the atom great 
compared with the wave length varies directly as the acceler- 
ation of the vibrator, which is proportional to its displacement, 
and inversely as the distance. The coefficient A is a function 
of the direction cosines #, /3, 7, of the line OP , as well as of 

*0’ V 

The total electric intensity at P, then, is equal to 




%Ae l 


i — ^ + + ct) 


E 




since r^E — (xa + y/3 + 27) 

is a sufficiently close approximation for the exponent. 

Put F= — {a l9 (<* 0 "- a )+ a iz(yo— 7)}’ 

A, 


H=^L {a nx (a^ -«) + a su (j3 0 - /3) 4- a 9z (y 0 - .7)} . 

A. 
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Then, by virtue of (18), 





+ n G + P H ) 

m n p 


(R-ct) 


It 




t ipII 


( 20 ) 


Hence the intensity of the diffracted radiation at P is maxi- 
mum for those directions for which 


X e ' mF ’ 2 /'"°’ X e ‘ pn 

m n p 

have their greatest absolute values. 

Suppose the illuminated portion of the crystal to be bounded 
by planes parallel to the sides of the elementary parallelepiped. 
Then m varies from — JR to + . 1 /, n from — N to + A, and p from 
— P to +P, where J If, N, P are positive integers. Hence 

X 4 mF = 1 + 4 F + e 2iF H b e mF 

* + c -«- + . .+e~ mF 

= 1+2 (cos 2'’+ cos 2F+ • • • + cos 1 WF}. 


The absolute value of this expression is obviously a maximum 
when F = ±2q 1 ir, 

where q l is a positive integer. 

Therefore the conditions for maximum intensity are 
a lx a + a ly ft + «,// = a lx a 0 + a tv ft 0 + a lz 7o ± (Jl \, 
a., x « + a iy ft + a . lz 7 = « 2x « 0 + a 2y ft 0 + a iz ± q 2 X, (21) 

a s x a ~b *3)//^ "b ~ a ix a o d~ a ayft 0 "b zV 0 ± 
where q and q 3 are also positive integers. 

The left-hand member of each of these equations is propor- 
tional to the cosine of the angle between one of the sides of 
the elementary parallelepiped and the diffracted ray, while the 
right-hand side is equal to an integral number of wave lengths 
plus a quantity proportional to the cosine of the angle between 
the same side of this parallelepiped and the incident ray. So 
for each value of q , q 2 , and q g these equations define three cones 
in a, ft, and 7. The loci of maximum intensity will be the 
lines of intersection of these cones. 
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Consider the simple case where a x , a 2 , and a 3 are mutually 
perpendicular and equal. Choose axes parallel to the sides of 
the elementary parallelepiped, and consider an incident wave 
advancing along the X axis. Then the con- 
ditions contained in (21) reduce to 
aa = a ± q t \, 

( 22 ) 

ar i = ± q 8 a. 

If a screen is placed at right angles to the 
X axis on the side of the crystal opposite 
to the source, the intensity of the trans- 
mitted radiation will be greatest along the 
traces on the screen of the cones defined 
by (22). The trace on the screen of the 
cone defined by the first of these equations 
is a circle, while those defined by the second and third are 
hyperbolae. The greatest intensity of all will be produced at 
those points on the screen where two, or better three, of these 
plane curves intersect. 

Consider radiation incident on the crystal under consideration 
at a small glancing angle 6 (Fig. 14). If the wave normal is 
in the XY plane, a o = sin0, 

/3 o = cos0, 

7 0 = °» 

and aa = a&in0 ± qfa 

a@ = a cos0 ± q 2 \, 

7 = ± ? 8 X - 

Consider the reflected radiation. For this 

a = — sin 6 , 
fi = cos dj 
7=0, 

and 2 a sin 6 = ± (23) 

? 2 = o, 



for reinforcement. 



CHAPTER VII 


ELECTROMAGNETIC FIELDS IN MATERIAL MEDIA 


37. Fundamental equations. Consider a medium composed 
positive and negative electrons. Experimental observation t 
the field in such a medium is limited, by the coarseness of avail- 
able instruments, to the investigation of the average values of 
the electric and magnetic intensities. These average values are 
defined in the following way. Divide the medium up into fixed 
elements of volume t so large that each contains very many 
electrons, but yet so small that no measuring instrument which 
may be used to investigate the field can detect a variation of 
electrical properties from one point to another in any one of 
these elements. Then the average value at any point of a 
scalar <f> which depends upon the state of the electrons in the 


medium is defined by 



where the integral is taken through the volume element t. 
Similarly, the average value of a vector V is defined by 



In the absence of an impressed field, each volume element of 
the medium will be supposed to contain equal numbers of posi- 
tive and negative electrons moving about in a fortuitous manner. 
Consequently the average charge and average current will 
vanish. In the presence of a field, however, the electric inten- 
sity may cause electrons of opposite sign to be displaced in 
opposite directions, with the result that the average density of 
charge may no longer be everywhere zero. In the same way, 
the magnetic intensity may orient intra-atomic rings of electrons 
in such a way as to produce an average current different from 

89 
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lalytical difficulties involved in discontinuities in the 
may be avoided by imagining every charged surface to be 
,d by a region of finite but very small thickness, in which 
es rapidly but continuously from one value to another, 
s the elements of volume r are fixed, 


OX TJ r OX 



etc., and similarly for the vector function V. 

To find the average value of p, consider a volume element t 
of dimensions Ax , A y, A z. In the time (it the charge entering r 
through those sides of this volume element which are perpen- 
dicular to the X axis is 

pv x AyAzdt - [pv x + £( pv x ) Azj AyAzdt 
= - faff®. r) Ax AyAzdt. 


Therefore, taking the six sides of t into account, 


Put 



— — V*(/»v). 



Then P = -V-Q, (1) 

and py = Q. (2) 

Hence, for a material medium, the equations of the electro- 
magnetic field given in section 26 take the following form : 

V-(E + Q)=0, (3) V-I = 0, (5) 

VxE=-*h, (4) V X H = ^ (E + Q), (6) 
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The method used in deriving these equations has not been 
such as to limit their applicability to media which are homo- 
geneous or isotropic. They apply to all media made up of 
electrons whatever their nature. Moreover, as two media in con- 
tact may be considered equivalent to a single non-homogeneous 
medium containing a thin transition layer in which the prop- 
erties of the medium change rapidly but continuously, 1 these 
equations may be applied in the region of contact. 

In order to determine the relation between Q on the one 
hand and E and H on the other, it is necessary to consider the 
motion of individual electrons in the medium. These electrons 
consist of two classes : (a) the free electrons, which move among 
the atoms, and (b) the bound electrons, whose displacements are 
limited by the boundaries of the atoms to which they belong. 

As an electron is of very small dimensions compared with an 
atom, the number of times a free electron collides with another 
electron is negligible compared with the number of times it col- 
lides with an atom. Hence if an electron strikes against an atom 
v times a second, its average drift velocity is given by 



where f is the average acceleration produced by the intensities 
E t and H 1 due to all charges other than that on the electron 
under consideration. Now 

mi = e^Ei-b - v X B^j 


approximately, and therefore if v is small compared to c the 
current due to the free electrons is given by 


— Ae 2 ■= 

py = o — E i 
2 mv 


where N is the number of free electrons per unit volume. 

Describe a sphere of volume 1 /N around the electron under 
consideration. Then if E a is the average electric intensity within 
this sphere due to this electron, 

E| “ E — E a . 
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In general E, is small compared to the other quantities 
involved in this equation. Hence (7) may be written 


pv = 


Ne 1 
2 mv 


1 


= C% (8) 

where C is the conductivity. 

A bound electron is supposed to be held in the atom to 
which it belongs by a force of restitution proportional to its 
displacement R from the center of the atom. In addition, it 
is supposed to be subject to a frictional resistance proportional 
to its velocity. Hence, if E 1 is the electric intensity and H x 
the magnetic intensity due to all causes external to the atom 
under consideration, the equation of motion of a bound electron 
inside this atom is 

?nR -|- aR -f- /SR = e d" ~ R X ’ 

where those terms in the dynamical reaction other than the one 
due to the mass have been omitted as negligible. 


Put 


2 Zee 


a 


m 



Then the equation of motion becomes 

R + 2ffi + i t 'R = -fE 1 +-RxH 1 l. 

m l o J 


(9) 


The displacement R satisfying this equation consists of the 
solution R 0 of the complementary equation 

R 0 +2lR t +k°R t =£R t xH 1 (10) 

added to a particular solution of the original equation. 

Obviously represents the part of the displacement pro- 
duced by the electric field E lt while R t is the part due to the 
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natural motions of the electrons inside the atom as modified by 
the presence of the magnetic field Since is in general 
negligible compared to c, the particular solution may be obtained 
to a sufficient degree of accuracy from the equation 

R,+ 2 JR, + (11) 

Now if there are N atoms per unit volume, and each of these 
contains n electrons, the current produced by the bound electrons 

is given by -7- -r— 

pv = NneR p + Nne R 0 . 

While eR Q may vary greatly from one electron to the next, 
it is evident that eR p will have more or less the direction ef 
the average electric intensity. Consequently the contribution to 



due to the entrance of new electrons into the volume element r 
through which the average is being taken will be negligible com- 
pared to the part of this quantity dependent upon the electrons 
already inside this region. Therefore 


and 


eR p = eR p , 

pv = NneR p + NneR Q . (12) 


If the electrons inside the atom are subject to different re- 
storing and resisting forces, this equation must be replaced by 


pv = (13) 

t t 

where denotes the number of electrons of type i inside 
each atom. 


The electric polarization, or electric moment per unit volume , is 


defined as 


P = Nne R p , 


or, if each atom contains more than one type of electron, 
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Hence (13) may be written 

pv^P+N^n^eR,. (14) 

At first sight it would seem as though eR 0 must always vanish. 
In a magnetic substance, however, this is not necessarily the 
case, as will now be shown. 

Suppose that for certain types of electrons the damping 
coefficient l in equation (10) is zero, and that the motion is 
constrained to a plane. Then motion in a circle with con- 
stant frequency is a solution of this equation. A substance 
each of whose atoms contains one or more rings of electrons 
revolving about their respective axes is said to be magnetic . 
The electrons in each ring will be assumed to be evenly spaced, 
and to have a constant angular velocity 0 about the axis 
of the ring so long as the external magnetic field remains 
unchanged. 

In the presence of a magnetic field each electron will be 
subject to a force which will tend, in general, to change the 
radius li Q of the ring and the frequency of revolution. More- 
over, the ring as a whole will be acted on by a torque which 
will tend to orient it. The force on an electron resulting from 
the magnetic intensity H x due to all causes external to the 
atom in which it lies is given by 

K = -(tx H ) 

C 

= ^ (fl x R 0 ) x H t , 

as R # is measured from the center of the ring. The torque about 
thi» point is g=2 R.XK 

for the entire ring. The form of this expression shows that G 
is perpendicular to 0. Denote by <f> the angle which R # makes 
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with the plane of 12 and H 1 . Then if 0 is the angle between 12 
and H 1? the component 6r„ of the torque parallel to the plane 
of these two vectors is 

G n — - sin 6 sin <j> cos <f> 


and the component G L at right angles to this plane is 
G l = ~ sin 0 cos 2 <£ 

= Sind, 

2 G 1 

where n is the number of electrons in the ring. The quantity 

2 c 

is called the magnetic moment of the ring. Therefore the torque 
is given by G = MxH, (15) 

Denote by dN the number of atoms per unit volume in which 
the axes of the rings of electrons have a given direction. Then 
the intensity of magnetization , or magnetic y 
moment per unit volume, is defined as I 


M dN, 


or, if each atom contains more than one Q P 

rin S’ I^yi, 

Consider a ring of electrons whose center is located on the 
X axis (Fig. 15) at a distance p to the right of the origin. Denote 
by 0 the angle which M makes with the X axis, and by <f> the 
angle which the MX plane makes with the XY plane. It is 

desired to find the value of • 

2> R o 
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due to electrons in this ring which lie to the right of the YZ 
coordinate plane. Evidently this sum will differ from zero only 
if the ring is cut by this plane ; that is, if 
— I£ o sin0 <jp <I2 o sin0. 


Moreover, it is obvious that this sum will have no component 
in the X direction. If yjr is the angle between R 0 and the inter- 
section QP of the plane of the ring with the MX plane, 


QH) v = sin cosi/r 


integrated from 
to 


Therefore 


( 2 *).- 


• ne B a 

2 7 r 

sin (f> J* cos 

— cos“ ] 

i( P \ 

\B. sin<?/ 

cos -1 

( P \ 

\R Q sm0/ 

nefl V Rl sin 2 # — p 2 

7 r 

sm 6 


Summing up over a prism of unit cross section extending from 
p = — E q sin 0 

to p= B o sin 

the y component of the current to the right of the YZ plane due to 
all rings which are cut by a unit area of this plane is found to be 


2 ( 2 *).=^ ■f^dfrJ^K;mr‘e- p ‘dp 

= 6 * J M&md sin <f)dN 




= cl. 


MAN 


Similarly, 


X(X e:k o\=-oi y , 


and, in general, 2(2*)— as) 

where n is a unit vector normal to the surface on the positive 
side of which the current is to be computed. 
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Consider a volume element t of dimensions Ax , Ay, A z. The 
current contained in it on account of those rings of electrons 
which are cut by the sides perpendicular to the X axis is 

I+^Aa^xiAyAz 

.OX / 

- < ’{- i S +k S} A * a! ' 4 * 

Therefore, taking the six sides of t into account, 

NneR 0 = cV X I, (17) 

or, if each atom contains more than one ring of electrons, 

iV^eR^cVxI. (18) 

t 

Hence the entire current due to bound electrons is given by 
pv = P + cVxL (19) 

Consider a medium in which there may exist currents due 
both to conduction by free electrons and to the displacement 
of bound electrons. Then it follows from (8) and (19) that 

Q = CE + P + cV x I, (20) 

and equation (6) becomes 

Vx(H-I) = i(i + P + C'E). (21) 

In accord with the usual practice put 

D = E + P, 

and B = H. 

The vector D is called the electric displacement , and B the 
magnetic induction. The average magnetic intensity due to all 
causes other than the intensity of magnetization I of the medium 
is usually denoted by H. As this letter has been used with 
another meaning in the preceding pages, this quantity will be 
designated by L. Then, omitting strokes, equation (6) gives 

V x L = - (E +P + CE), 

c 

and comparison with (21) shows that 

L = B — I, 
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since L vanishes with B and I. As B is the total average mag- 
netic intensity, I represents the average magnetic intensity 
produced by the magnetization of the medium. Equation (21) 
may now be written 

V x L = 1 (f> + CE). (22) 

c 

Returning to (20), it follows that 

V*Q = CV*E+V*P, 

and V-Q = C fv-Edt+V-P, (23) 


as Q and P vanish everywhere when no external field is present. 
Therefore, omitting strokes, equation (3) becomes 


• d ~ c / 


V«E dt. 


Hence the equations of the electromagnetic field in a medium 
containing both free and bound electrons take the form 


V-D =~cj V-E dt, (25) V-B = 0, 


VxE = 


VxL = ^-(D + (7E), (28) 


when the field is investigated from the macroscopic point of view. 
As already noted, these equations apply whether or not the 
medium is homogeneous or isotropic, and they are valid in the 
region of contact of two media. 

If, in addition to the processes of conduction and displace- 
ment, charges p and currents pv are produced by convection 
through the medium, the electromagnetic equations for this most 
general case assume the form 


V'D = ~cJV'Edt+p, (29) V*B = 0, 


V x E = — - B, 

c 


(30) VxL=i(D+C , E+pv). (32) 
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38. Specific inductive capacity. In order to obtain the relation 
between P and E, and hence between D and E, it is necessary to 
solve equation (11) of the preceding section. The following dis- 
cussion will be confined to the case of simple harmonic fields, 
a steady field being considered as a simple harmonic field of 
zero frequency. The electric and magnetic intensities are given 
by the real parts of the following complex quantities, 

E = E 0 <r tW , 

, d 

and ~r = — ico. 

dt 

Hence the solution of (11) is 

e 

T* m T> 

■* !=S yfc 0 2 — cw 2 — 2 iQ)l El ’ 

Nn — 

7YI —— 

and P is given by P = E r (33) 

To find the relation between E A and E, describe about the 
center 0 (Fig. 16) of each atom a sphere of volume 1/K Then 
if E 2 is the average electric intensity 
within one of these spheres due to the 
electrons which it contains, 

E x = E — E 2 . 

Evidently E a has the direction of 
the displacement caused by the ex- 
ternal field Ej, and its magnitude is 
given approximately by Fig. 16 

E 2 = NnCC " - - 2 cos 6 (2 i rr 2 sin OdOdr ), 
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where 0 is the angle which the radius vector makes with R^. 
As r goes from zero to 



— COS 0 + COS 2 # + ~ — 1 1 » 


it follows that 

E 2 = — ^ NneRp 

(34) 

Therefore 

E^E + iP, 



Nn — 


and (33) becomes 

r» m _ 

lc 2 — © 2 — 2 icol ’ 

(35) 

where 

tfEzk*-\Nn-. 

3 m 


Put 

Nn — 
m 

k 2 — co 2 — 2 zW 



Then, omitting strokes, 


P = eE, (86) 

and D = (l+e)E. (37) 

The specific inductive capacity k is defined as the ratio of the 
displacement D to the electric intensity E. Hence 


k = 1+ e 


= 1 - 


Nn — 
m 


ti l — co 2 — 2 icol 


(38) 


or, if each atom contains more than one type of electron, 

* sl +S e< 


= 1+ N— V ^ — 7 

m k 2 — co 2 — 2 ml. 

In the case of a steady field, 


k = 1 + 


Nne 2 
k 2 m 


(39) 


(40) 


if there is only one kind of electron in each atom, and 
Nhe 2 fc — 1 

3 k'£m fc + 2 


( 41 ) 
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The number of atoms per unit volume varies directly with 
the density d. Therefore ^ ^ 

<i(# + 2) 

should be a constant. This relation has been verified for many 
gases. 

The force of restitution in a polarized atom may be easily 
evaluated if n positive and n negative electrons are assumed 
to be distributed uniformly over equal and coincident spheres. 
If the center of one of these spheres remains stationary, while 
that of the other is displaced a distance the force drawing 
them together is approximately 




where v is the volume of a sphere. Hence the force on each 
electron is .. 2 


and 


1 ne 2 
3 vm ’ 


Substituting in (41), it is found that 


Nv = 


K — 1 

tc + 2 


(43) 


is the portion of the volume under consideration which is 
actually occupied by the atoms. If an estimate of JV has been 
obtained from other sources, the atomic volume v may be com- 
puted. In this way the radius of an atom is found to be of 


the order of 


(10)“ 8 cm. 


39. Magnetic permeability. It has been pointed out in section 37 
that a magnetic field tends both to change the radius and the 
frequency of revolution of the ring of electrons inside a magnetic 
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atom, and to orient this ring. To determine the magnitude of 
the first effect, integrate both sides of the equation 

V x E = — - Hjt 
e 

over the surface bounded by a ring of electrons. Then 

Vx E^da = — -J* H ^da. 

But by Stokes’ theorem 

J * V x E 1 *c?<r = J* E^c^X, 

where the right-hand member is integrated around the ring. 
Hence the work done on the ring when the magnetic field is 
increased by dH is n ~ 

It I E fd\ 


d(J = ™*} L dt 

A IT 


= _!£^ a .H M 
2c 1 


(44) 


Now the increase in the kinetic and potential energies of the 
electrons in the ring is 

dU = nm(ydv + k^R^dR^) 

= nmRjCl ( R 0 dH + 2 D<£K 0 ), 


as 


But 

Therefore 


k 0 — n. 

dM= ^ ( R 0 <m + 2 

A C 

, rr 2 mSle J JT 
dU = dM. 


(45) 


Equating this to the expression for the rate at which work 
is done on the ring, 

dM=- — M- dR v (46) 
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Let a be the angle which M makes with c?H 1 . If initially the 
axes of as many atomic rings pointed in one direction as in any 
other, dl has the direction of dH 1 and is given in magnitude by 



dMao&adN 


eMdU l r 
2 mtlcj 


cos 2 adN 


NeM 
6 m£lc 


dH t . 


(47) 


To find the relation between dTl x and dH, describe about the 
center of each atom a sphere of volume 1/N. Then if H a is the 
average magnetic intensity within one of these spheres due to 
the electrons which it contains, 

H, = H — H 2 . 

But it has been shown that the average magnetic intensity H 2 
due to rings of electrons is equal to the intensity of magneti- 
zation I. Therefore as b ^ H 

it follows that H x = B — I 

= L, 

and dl = --^-dL. (48) 

b mile 


The permeability p of the medium is defined as the ratio of 
the magnetic induction B to the external field strength L . 
Therefore in the case under consideration 


NeM 
6 m£lc 


(49) 


The effect under discussion is known as diamagnetism , and is 
characterized by a value of the permeability less than unity. 
It is shown in the greatest degree by bismuth. 

The tendency of a magnetic field to orient a ring of electrons 
will be considered next. The torque on such a ring has already 
been shown to be n wTT 
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Let H be the average total magnetic intensity inside a sphere 
of volume 1 /N described about the center of the atom under 
consideration, and let H 2 be the average magnetic intensity 
within this sphere due to the electrons which it contains. Then 

h^h-h,, 

where H 2 evidently has the direction of M. Therefore 

G = M x H. (50) 

Consider a ring of electrons which is rotated from a position 
where M is perpendicular to H to one where M makes an angle a 
with H. The potential energy acquired is 

U=MII Ji sin ada 

= — MU cos a 

= — M«H. (51) 

Therefore if the tendency of the magnetic field to bring all 
the atomic rings into line is opposed only by the disorganizing 
effect of thermal agitation, the number of rings per unit volume 
whose axes make angles with H between a and a 4- da is given by 

u 

dN~\Ae * kT sin ada 

MH cos a 

= \Ae * kT sin ada , 

where T is the absolute temperature and \hT the average 
kinetic energy associated with each degree of freedom. 
Integrating, the constant A is found to be given by 

Nx 

sinh x 3 


where 


MU 



MB 

| kT m 
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The magnetic moment per unit volume has the direction of B 
and is given in magnitude by 



MvosadN 



cosh# 
sinh x 



(52) 


As B increases I does also, approaching the saturation value 
NM for large values of the field. For small fields 


/ = ~ NMx 
o 


2 NM 2 
9 kT 


B, 


(53) 


showing that the intensity of magnetization is proportional to 
the strength of the field. The permeability is given by 


, , 2 NM 2 
+ 9 kT ' 


(54) 


The effect under discussion is known as paramagnetism, and 
is characterized by a value of the permeability greater than 
unity. It is shown to an exceptional degree by iron, in which 
the effect is given the special name of ferromagnetism. The theory 
given above cannot be considered as more than very roughly 
approximate to the facts, especially as it gives no explana- 
tion of hysteresis, or the lagging of the magnetization behind 
the field. 

40. Energy relations. From (30) and (32) it follows that 

L*VxE-E-VxL=--(E*D + L*B) — - (CE 2 + pE*v). 

c c 

But L«Vx E-E»VxL = V*(Ex L). 

Hence 

j t {- (kE 2 + pZ 2 )} + «V- (ExL)+ CE 2 + pE-v = 0. 
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Integrating over any arbitrarily chosen portion of space, and 
applying Gauss’ theorem to the second term, 

J t j\ (*-® 2 +/*£*) dT + c J"(ExL) -d<r+ (CJ^+pE-v) dr=0. (55) 


The third term of this expression measures the rate at which 
work is done by the electromagnetic held on the conduction 
and convection currents in the medium. Following the line of 
reasoning pursued in section 28, the conclusion is reached that 
the first term represents the rate of increase of energy of the 
held, and the second the hux of energy through the surface 
enveloping the held. The forms of these expressions suggest that 

u = 1(k!£*+iiL*) 


is to be considered as the electromagnetic energy per unit 
volume, and 8 = c(E xL) 


as the how of energy per unit cross section per unit time. 

41. Metallic conductivity. In developing the electron theory 
of metallic conduction, the atoms in a metal may be treated as 
immobile compared with the electrons. Conduction currents of 
electricity in a metal will be supposed to be due entirely to the 
drift velocity of the free electrons in the direction of an im- 
pressed electric held, and heat conduction will be attributed to 
the transport of energy by these electrons from one atom to 
another in the direction of the temperature gradient. 

If u is the average velocity of an electron due to thermal 
agitation, and l the average path described between collisions 
with atoms, the number of times an electron strikes an atom 
per second is given by 



and equation (8) for the conduction current becomes 


C E = 


M. B. 

2 mu 
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As | mu 2 = | kl\ 

the electrical conductivity is given by 


NeHu 
6 kT ' 


(56) 


In determining the heat conductivity, take the X axis in the 
direction of the temperature gradient. Consider an electron 
which is just about to collide with an atom at a distance x from 
the origin. This atom will have, on the average, a kinetic energy 

but the electron will have the energy 

Icond) 

2 \ dx ) 

of the last atom with which it. collided, where 0 is the angle 
which the electron’s path makes with the X axis. During the 
collision the electron will come into thermal equilibrium with 
the atom, giving to the latter an amount of energy equal to 

S 7 dT 7 a 
2 dx 


Now the number of electrons per unit volume whose paths 
make angles between 6 and 0 + d6 with the X axis is 

|iVsi n0d0, 

and the number of these which pass through unit area at right 
angles to the X axis in unit time is 

1 Xu cos 0 sin 0d0. 


Hence the flux of energy is 


dT 


-Nluk — 
4 dx 


* r*n 
1 ' 


cos 2 0 sin 0d0 


'muk 8 * 

2 dx 


giving for the thermal conductivity 

K= i Xluk. 


( 57 ) 
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The ratio of thermal to electrical conductivity, 



( 58 ) 


varies directly with the absolute temperature and is the same 
for all metals at a given temperature. Its value for any tem- 
perature depends only upon the universal constants k and e , 
whose values may be determined from experiments having no 
connection with the metallic conductor in question. The ratio 
of the conductivities as thus computed is in fair agreement 
with the ratio as determined directly by experiment. 

In obtaining the expressions for the conductivities given 
above, use has been made of the average velocity of thermal 
agitation and the average length of path between successive 
collisions. A more exact calculation gives a slightly different 
numerical coefficient for the ratio, but one which shows rather 
worse agreement with the experimental value of this quantity. 

42. Reduction of the equations to engineering form. If free 
charges and currents — either conduction or convection — are 
present in a material medium, the equations (29) to (32), 
section 37, of the electromagnetic field may be written in 
the form 

V*D = —J' V • J dt, (59) V*B = 0, (61) 

VxE = -*B, (60) Vx L=4(f> + j), (62) 


where J is the current density; that is, the current per 
section. Moreover, J = C E -f- pv 

D = E 4- P 


( 68 ) 


= *E, (64) 

B = L + I 

= fiL. (65) 

The quantities involved in these equations are measured in 

Heayiside-Lorentz units. If the same quantities as measured 
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in electromagnetic units are designated by letters with the sub- 
script m, and as measured in electrostatic units by letters with 
the subscript a, then 

p = c V4 irp m = V4 irp s , 

J = eV4wJ m =V4^J„ 

P = c V4 7rP m = V4 7rP s , 

I=V4 Trim — 6 *V4 7rl s , 

D = <?\/4 7rD m = V 4 ttD., 


B = — — B = 
/l — 171 

V 4 7T 


V4^’ 


E = 


V4 7T 
1 


E =s 


L = — — L = 
VTtt " 

C=\ir#C = 


1 

V4 7 r 
1 


E. 


V4 7T 

4 irC„. 


Therefore, in electromagnetic units the field equations take 
the form 

V*D m = -jV JJt, (66) V-B m =0, (68) 

VxE m =-B m , (67) VxL ni =4w(D m +JJ, (69) 

where J m = C m E m + p M v, (70) 


4 7T6* 5 


a E m + P„ 


K 

4^? 


E 




(71) 


= + 4 7rl m 

= /*L.- (72) 


In electrostatic units these equations are 
V-D ,= rf, (73) V.B i= 0, (75) 

VxE,= -B„ (74) VxL 5 =47r(i> 4 +J,), (76) 
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where the charge accumulated at any point has been denoted by 
p' instead of r 

-J v-j ,dt, 

and J s = C a E,s+ p,v, (77) 

D s =^-E s +P s 

: 7 T 


- K v 

“ 4 7T Es ’ 


(78) 


B,= ? L s +4wI s 


_ M t 

_ 7 l *' 


(79) 


In practical applications of electrodynamics these equations 
are generally made use of in integral form. For instance, con- 
sider a small charge e permanently at rest at the center of a 
sphere of radius r. Integrating (73) throughout the region 
enclosed by the spherical surface, 


But 


/ 7 

/ v 


•D s dr = 


-D a dr = 


f"' A 

e s* 

J B s -d<r 


= 4 1 rr 2 D, 

by Gauss’ theorem. Hence 


and, by virtue of (78), 


D , 


4 7 rr 2 


E> 


( 80 ) 


Again, integrating (67) over a surface, 
Jv X E m -efcr = -J B m ‘dor. 
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But, by Stokes’ theorem, 

JV* E m 'da = J E m -d\, 

where the line integral is taken along the boundary of the 
surface. Therefore # 

JE m -d\ = -jB m -dcr, (81) 

which is the usual form of Faraday’s law of current induction. 
Similarly, (69) leads to 

J L.-dX = 4 7r J ‘(ft.+jJ-dr, (82) 

which is the form in which Ampere’s law is generally expressed. 



CHAPTER VIII 


ELECTROMAGNETIC WAVES IN MATERIAL MEDIA 

43. Isotropic non-conducting media. For wave lengths long 
compared with atomic dimensions, the electromagnetic field is 
specified by equations (25) to (28) of the last chapter. For 
non-conducting media these take the form 

V.D = 0, . (1) V*B = 0, (3) 

V x E = — i B, (2) V x L = - D, (4) 

G C 

where D = vtE, 


Eliminating B as 

B = /jlL. 

in section 29, it is found that 



V.VE — *?E = 0, 

c 

( 5 ) 

or, eliminating D, 

V-VL - ^ L = 0. 

c 2 

( 6 ) 

These are equations of a wave travelling with velocity 



II 

*v 

( 7 ) 


The permeability of all substances is practically unity for 
frequencies as great as that of light. Hence for light waves it 
is permissible to write c 

v r= ‘ (8) 

V* 

The index of refraction n of a medium is defined as the ratio 
of the velocity of light in vacuo to that in the medium. 

Therefore w 2 = /c, (9) 

and, as was shown in section 38, 
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(10) 
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In the case of a plane wave it is often convenient to make 
use of the wave slowness S in place of the wave velocity V. 
This quantity is defined as a vector having the direction of the 
wave velocity but equal in magnitude to its reciprocal. Hence 
the index of refraction may be defined as the ratio of the wave 
slowness S in the medium to the wave slowness S 0 in vacuo. 
The electric intensity in the case of a simple harmonic plane 
wave advancing in a direction making angles a , /£, 7 with the 
X, F, Z axes, may be expressed in terms of the wave slowness 


by the real part of 


i£ = j0/ w(S - r -‘>, 


which is far more compact than the equivalent expression 

( a? cos + jy cob fl + z cos y \ 

y t) 

involving the wave velocity. 

Equation (8) shows that the wave slowness is given by 

S=S 0 ^K, ( 11 ) 

where k is in general complex. To show the significance of a 
complex wave slowness, put 

SsS'+tS", 

where S' and S " are real. Then 

E = E 0 e-” s "- r^S'-r-o (12) 

showing that the imaginary part of the wave slowness measures 
the damping of the wave as it progresses into the medium, 
whereas the real part determines the actual velocity of propa- 
gation. The same statement applies to a complex index of 

refraction. If , . N 

n==v(l + ix ), 

where v and % are real, it follows from (10) that 

(18) 


( 14 ) 
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For the portion of the spectrum to which these expressions 
are to be applied it will be assumed that 

*,■»<•',»*•. j- 
Then x^l, 

mi ?=l (15) 

approximately. Consider the denominator 
(k? — &> 2 ) 2 -f 4 o> 2 Z 2 

of one of the terms in the right-hand member of this equation. 
Except for the region where co is very nearly equal to k { the 
second term of this denominator is negligible compared to the first. 
Within this region the first is negligible compared to the second. 
So if v 2 is plotted against co 2 in a region extending from go 2 to co 2 , 
such that k x lies between aq and eo 2 , and fc 2 , /c 3 , etc. lie outside 
this range, 2 2 

Nn, - Nn~ 

■'- 1 + 5 = 5 + )?= 5 + *- (16) 

except in the neighborhood of h*, where 


Nn — ( k 2 — co 2 ) Nn — 


; -b etc. 


Plotting each term separately, the dotted curves of Fig. 17 
are obtained. Adding these curves, the full line curve is found 
to give the relation between v 2 and co 2 . The portion of the curve 
from A to B corresponds to regular dispersion, the index of 
refraction increasing as the frequency becomes greater, whereas 
the part BC accounts for the anomalous dispersion observed in 
the neighborhood of an absorption band. It must be remem- 
bered that the index of refraction refers always to the phase 
velocity of a train of waves. Hence the fact that this index 
becomes less than unity on the short wave length side of an 
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absorption band may not be adduced to show that an electro- 
magnetic signal can be despatched with a velocity greater than 
that of light in vacuo . In fact Sommerfeld has shown that the 



Fig. 17 


forerunners of a limited train of waves travel with the same 
velocity c in all media and are undeviated when they pass from 
one medium into another. 

44. Anisotropic non-conducting media. For wave lengths long 
compared with atomic dimensions, the first four equations of the 
electromagnetic field are the same as for isotropic media. The 
relation between D and E, however, is different, as the atoms in a 
body which is not isotropic must be supposed to exert different 
restoring forces in different directions. Hence equation (36), 
section 38, for the polarization must be replaced by the more 
general relation = ^ + e ^ a + e 13 E„ 

and similar expressions for P y and P z . In the following dis- 
cussion damping will be assumed negligible. Therefore the 
coefficients of the components of E will be real, though func- 
tions of the frequency of the radiation traversing the medium. 
In vector notation p _ (18) 

where is the dyadic eU+eij+eik. 

+ + € J k 

+ € 8 l ki + e 32 k j+ e 8 3 kk ' 
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If u denotes the energy per unit volume of the medium, 
du = E»dP 

= («U*. + e n E V + € 81^) dE x 

+ ( e n E x + e 22 E y + € 32 E z) dE V 
+ ( e n E * + e n E v + € 88 E z) dE z- 

Now the law of conservation of energy requires that this 
expression shall be an exact differential. Therefore 

• € ij 

and the dyadic \|f is self -con jugate. Hence by a proper choice 
of axes \| i may be put in the form 

^Esejl + ejj + ejtk. 

Now D = E + P 

= (1 + \|/*) E. 

Put $ ~ k x ii + K y ] j + ^kk, 

where ^==1+6*, 

K v = 1 + e z/> 

** = 1 + 6 *- 

Then D = $-E, (19) 

showing that the specific inductive capacity is a dyadic instead 
of a scalar factor as in the case of an isotropic medium. For 
D and E are not in general in the same direction in a medium 
which is not isotropic. Eliminating B and D from equations (1), 
(2), (3), (4), and (19), it is found that 

V*VE — VV*E — = 0. (20) 

Confining attention to plane waves, the electric intensity and 
displacement are given by the real parts of 

E = E/"(S- r -0, 

and D = D # e <w ( s,r_, >, 

respectively. 

Therefore d E = dr* (i»S) E — dticoE. 

dE = dfVE + dtE. 


But 
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Comparing, V = icoS, 

and ~ = — ico. 

dt 

Hence (20) becomes 

Sfi* E + SS-E - S*SE = 0. (21) 


Multiplying by S*, it is seen that 
S»D = 0, 

or the vector D is at right angles to the direction of propa- 
gation of the wave. Moreover, equation (21) shows that D, E, 
and S are in the same plane. 

The magnetic intensity may be found from (2). This equation 


gives 


S x E - N 0 L, 


showing that L is perpendicular to both S and E, and hence 
to D. The vectors D and L lie in the wave front, although E 
makes an angle with this plane. The flux of energy is given by 


s = cExL, 


which is a vector at right angles to E in the plane of D, E, and S. 
Fig. 18 shows the relative directions of .the vectors under 
discussion. Since the flux of energy is not along the wave 
normal, limited wave fronts will side- 
step as they advance, as indicated in the 
figure. A line drawn in the direction 
of the flux of energy is known as a ray, 
and the velocity of propagation of energy 
along this line as the ray velocity . 

Consider an infinite number of plane 
waves passing through the origin atjthe 
time 0 in all directions, the vector D 
having all directions in the wave front. 

The envelope of these plane waves one second later is known 
as the Fresnel wave surface. To find the equation of this sur- 
face, it is necessary to obtain from (21) a relation involving 
the wave slowness as the only unknown quantity ; that is, an 
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equation between S, S 0 , and $ which is true for all possible 
directions of E. Equation (21) may be written in the form 

(£ 2 $ + SS - S-SI) *E 0 = 0, (22) 

where I is the idemfactor 

ii + jj + kk. 


The dyadic in the parentheses in (22) causes the vector E 0 
to vanish. Hence either its antecedents or consequents must be 
coplanar. This dyadic may be written in the expanded form 

- $?) i + W + S.SJl} i 

+ {SJS, i + {S 0 \ - - ^ 2 ) j + s 9 s; k} j 

+ {s x s z i + s y sj + s?) k}k. 


As the consequents are not coplanar, the antecedents must 
be. Therefore the scalar triple product of the latter must vanish. 

Hence (sx -s 2 + s?) CSX -.s i 2 + sf) (s 0 \ -s*+ sfi 

- s;s: asx -s 2 + - s 2 s 2 (s 0 \ -s 2 + 

- S?s* (SX -s 2 + S?) + 2 s;s;s; = o, 

or, reducing. 


1 + 


S? 


+ 


s? 


s? 


sx-s 2 sx-s 2 sx-s 2 ' 

If l, m, n are the direction cosines of the wave normal, this 
equation becomes 


l 2 


m 


n 


V 2 ~a 2 V 2 -V V 2 -<? 


0 , 


( 28 ) 


where V is the wave velocity, and 
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The equation of a plane wave advancing in the direction 
specified by l, m, n is h + m;/ + m _ r (!4) 

one second after leaving the origin, where Z, wi, n satisfy the 
condition P+™*+„'_l. (25) 


The Fresnel wave surface is the envelope of the family of 
planes obtained by varying Z, m , n in equation (24), subject to 
the conditions specified in (23) and (25). This surface is most 
easily found by differentiating (23), (24), and (25), and elimi- 
nating Z, m, n, and V by means of the original equations. The 
equations obtained by differentiating are 



xdl 4- ydm + zdn = dV, 

(26) 


Idl + mdm 4- ndn = 0, 

(27) 


idl mdm ndn , , 

F 2 — a 4 1 F 2 — b 2 ' V 2 -c*~ ' 

(28) 

where 

j _ l 2 to 2 n 2 

(F 2 -a 2 ) 2 1 (F 2 — 6 2 ) 2 (F 2 — c 2 ) 2 * 

(29) 


Let — p and — q be factors by which (27) and (28) respectively 
may be multiplied so as to eliminate dl and dm when the three 
equations are added together. As the other differentials are 
independent of one another, the coefficient of each differential 
in the sum must vanish, giving 

X =P l 


m 

y=pm + qys-j 2 , 

z =P n + 5 ^ 37 ’ 

i=qkr. 


(30) 


(31) 


Multiplying equations (30) by l, m, n respectively, adding, 

and making use of the original relations (23), (24), (25), it is 

found that T . , 00 , 

p — V. (32) 
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Squaring and adding the two sides of equations (BO), 

r* = p i +q% 

the product term disappearing on account of (23). Combining 
with (31) and (32) q = (f-V^V. (33) 


Substituting in (30) the values of p and q just found, 

V(t>-V)l VI (r* — 
x — ' * ' V 2 cP V 2 cv 2 * 


X VI x-Vl 

r*-a 2 ~~ V 2 -a 2 ~~ t*-V 2 ' 
V __ V' m y — Vm 

T^b 2 ~~ V^-b 2 ~ r 2 — V 2 ’ 
z _ Fn __ 2 — Ftt 
r 2 -c 2 ~~ ~V 2 — <? ~~ r 2 — V 2 * 



Fig. 19 



Multiplying these equations by 2 respectively and adding, the 

required equation ^2 ^2 ^ ^ ^ 


+ 


+ 


r 2 -b 2 ' r 2 — c 2 
of the Fresnel wave surface is obtained. 

In discussing this surface, it will be considered that 

a 2 <b 2 < c 2 . 

The trace of the surface on the YZ coordinate plane is 
expressed by 




which is a circle inside an ellipse, as in Fig. 19(a). On the ZX 
plane the trace is /?* ji \ 
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which is a circle cutting an ellipse, as in Fig. 19(6). Finally, 
the trace on the XY plane is given by 



which is an ellipse inside a circle, as in Fig. 19 (<?). 

One octant of the surface is represented in Fig. 20 ( a ), and 
a section through the point P is shown in 20 (6). As D lies in 
the plane of the wave normal and the ray, this vector must be 



tangent to each of the elliptical traces, as indicated by arrows 
in the combined figure. Consequently it must be perpendicular 
to the planes of the circular traces. 

The primary optic axes of a crystal are defined as those direc- 
tions in which the wave velocity is independent of the state of 
polarization, that is, the direction of D in the wave front. Hence 
the perpendicular OQ to the tangent QT [Fig. 20 (a)] is one 
of the primary axes. The other is also in the ZX plane, making 
an equal angle with OX on the other side of the X axis. A 
crystal which has two optic axes is known as biaxial. Obviously, 
there can be no more than two such axes. Since the constants 
which determine the intercepts of the Fresnel surface depend 
upon the three principal specific inductive capacities, which are 
themselves functions of the wave length, the directions of the 
optic axes of a biaxial crystal vary with the wave length. 
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The secondary optic axes are defined as those directions in 
which the ray velocity is independent of the state of polari- 
zation. One secondary axis has the direction OP , the other 
making an equal angle with OX on the other side of the X axis. 

A uniaxial crystal is one in which two of the quantities #, 6, c 
are equal. If b and c are equal, the crystal is said to be positive 
or prolate . The Fresnel wave surface is shown in Fig. 21 (a). 



There is only one axis, and there is no longer any distinction 
between primary and secondary axes. Moreover, the direction 
of this axis is independent of the wave length. 

If a and l are equal, the crystal is said to be negative or 
oblate . The Fresnel wave surface for such a crystal is shown in 
Fig. 21 (5). 

45. Reflection and refraction. Consider a train of plane waves 
incident at an angle' (Fig. 22) on a plane surface separating 
two transparent isotropic media. The incident light will be 
partly reflected and partly transmitted. Let A x be the amplitude 
of the electric vector in the incident radiation, A[ that in the 
reflected, and A 2 that in the transmitted radiation. Then the 
coefficient of reflection R is defined by 



and the coefficient of transmission T by 
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Let the subscripts 1 and 2 refer respectively to the media 
above and below the plane 0 Y. In the case of the upper medium, 
letters without primes will refer to the inc?ident light, and letters 
with primes to the reflected light. Attention will be confined 
to electromagnetic radiation of wave length long compared to 
the distances between adjacent molecules of either medium. 
Hence equations (1) to (4) inclusive specify the field. Moreover, 
if the media under consideration are transparent, the damping 
term in expression (10) for the specific 
inductive capacity is negligible, and 
this quantity is real. 

Choose axes as indicated in the 
figure, the Z axis extending upward 
from the plane of the paper. Consider 
a short pill-box shaped surface, with 
bases parallel to the YZ plane and axis 
bisected by this plane. Integrating 
(1) over the volume enclosed by this 
surface, and transforming the volume 
integral into a surface integral by 
means of Gauss’ theorem, it is found 
that the components of D normal to the surface of separation 
are the same on the two sides of this surface. A similar relation 
between the normal components of B follows from (3). 

Consider a rectangle of which one pair of sides is very much 
longer than the other, so situated that the short sides are perpen- 
dicularly bisected by the YZ plane. Integrating (2) over the 
surface bounded by .this rectangle, and transforming the left- 
hand side of the equation into a line integral by Stokes’ theo- 
rem, it is found that the components of E parallel to the surface 
of separation are the same on the two sides of this surface, 
provided B is not infinite at the surface. A similar relation 
between the parallel components of L follows from (4). 

Suppose the electric vector in the incident wave to be perpen- 
dicular to the plane of incidence ; that is, the light is polarized 
in the plane of incidence. Then the x and y components of the 


l 
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electric intensity are zero for each wave, and the z components 
are given by the real parts of 

glut { S t (x cob <J> i + 

I f> iu > { s i (- ^ cos t/),+ ?/ Bin </>,)- 1 

E = k A { S » ( XC0B ^2 + " ein 4>i) - *} 

2 2 “ * 


Hence, remembering that n~ 



= cS, 


and that the permeability is unity for the frequencies under 
consideration, it follows from (2) that 


L 1 = iA’ 1 n 1 sin cf> 1 — j E 1 n 1 cos 
LI = i E[n x sin <f> x -f j E[n x cos 
L s = sin ^ cos <f> 2 , 

showing that L x = E x n v L[ = E[n v and L 2 =E 2 n 2 . 

Therefore the relations between E and L on the two sides of 
the surface of separation lead to the three equations 

E 1+ E[ = E^ ■ ( 85 ) 

0 E , + K) «, sin </> I = A> 2 sin 0 2 , (36) 

i E i ~ E [) n i cos <#»i = ■»,», cos <f> 2 . (37) 

From (35) and (36) it follows that 


sin 

s>n \ S t ’ 


(38) 


the familiar relation attributed to Snell. Dividing (37) by (36), 


E x — E[ __ sin <£, cos </> 2 
E x + E[ sin <f> 2 cos cj) x 


(39) 


Now the exponentials in E t and E[ are the same when x is 
zero. Hence the electric intensity may be replaced by its ampli- 
tude in this equation. Thus 


1 — JKjl __ sin <f> t cos cf> 2 
1+jK ± sin^cos^ 
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determines the coefficient of reflection Jt ± for the case where the 
electric vector is perpendicular to the plane of incidence. Solving, 

B = shl ^ CQS ^ ~~ sin ^ CQS ^ = sin ~ M . MO) 

“ L sin $ 2 cos <f> 1 -f sin <^> 1 cos <f> 2 sin + <t> 2 ) * ^ ^ 

To find the coefficient of transmission T ± , eliminate E[ between 
(35) and (39). Since it follows from (38) that the exponentials 
in E x and E 2 are the same when x is zero, the electric intensity 
may be replaced by its amplitude, as before. Hence 
2 — T ± __ sin <j> i cos <f> 2 
1\ sin (f> 2 cos </> 1 ’ 
m 2 sin <f> cos 6 , 

whence T , = — — — - 1 . (41) 

In the case where the electric vector is in the plane of incidence, 
E 1 = {— i A t sin c{> 1 + ]A 1 cos , 

El= {— i A[ Sin </> 1 — j A[ COS gMM-aco.fc + ysM,)-*}, 

E 2 = {- i A 2 sin <f> 2 + jA 2 cos <f> 2 } e iul 1 
and it follows from (2) that 

L 1 = kA> 1 , 

L;=k^\, 

L t = kE 2 n 2 . 

In this case the relations between E and L on the two sides 
of the surface of separation lead to the three equations 

sin </>,=-&>* sin <j> 2 , (42) 

(2i\ - Ji[) cos cf> i = E 2 cos (48) 

(l^ + E'Jn-E^ * (44) 


These equations lead to Snell’s law and the following ex- 
pressions for the coefficient of reflection i£„ and coefficient of 
transmission T p 

jt = sin ^ cos ~ sin $■> cos & _ tan (jx ~ &} (45) 

11 sin <j> 1 cos <f> 1 + sin cos (f> 2 tan (<f> t + <£ s ) ’ 

T 2 sin 0, cos <£, 

" sin (<#>! + </> 2 ) cos (^>j — <£ 2 ) ' 


( 46 ) 



126 AN INTRODUCTION TO ELECTRODYNAMICS 


Examination of the four coefficients of reflection and trans- 
mission shows that R H is the only one which can vanish. The 
polarizing angle & x is defined as the angle of incidence for which 
this coefficient becomes zero. Therefore 

tan (3> 1 + <& 2 ) = oo , 

and <t> x and 4> 2 are complementary. Consequently 

Consider unpolarized light striking a surface at the polarizing 
angle. The energy associated with the component of the electric 
vector in the plane of incidence will be entirely transmitted. 
Consequently the reflected light will consist altogether of radia- 
tion in which the electric vector is at right angles to the plane 
of incidence. Although polarization by reflection should be 
complete at the polarizing angle, experiment does not show it 
to be so. This is largely due to imperfect surface conditions. 

If n x is greater than the incident radiation may be totally 
reflected. For convenience the discussion will be restricted to 

the case for which . 

n l = n > 1, 

n 2 = l. 

Then, for total reflection, 

and sin (j> 0 = n sin <j> 1 

is greater than unity, and 

cos <f) 2 =iy/ n 2 sin 2 ^ — 1 

is imaginary. Reference to (40) and (45) shows that both 
coefficients of reflection will be of the form 


a + ib 

E[ = A x Re ita { s i (° + v 8in *i) - *} 


Therefore 
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showing that total reflection has produced an alteration in phase 
without change in amplitude. 

Consider the transmitted ray in the case of total reflection. 

JU __ e iu { »?2 ( x cos </> a + ?/ sin <#>.,)-«} 9 

and as cos<f> 2 is imaginary, 

= A 2 e~ b} v x e i< ° IV Bin< fo- 
where iy = S 2 cos c/> 2 . 


This is a wave travelling along the surface, with an amplitude 
which falls off exponentially with x . As the wave does not 
pass across the surface, no energy is taken from the incident 
radiation. 

46. Rotation of the plane of polarization. Consider a beam of 
monochromatic plane polarized light travelling through a trans- 
parent medium in the direction of an impressed magnetic 
field H. Choose axes so that the X axis is parallel to H. Fol- 
lowing the method developed in section 38, it is found that 
equation (9), section 37, leads to 

e 

_ m f— .a) - 

R =F^n E- *c RxH r 


where the term involving the damping constant l in the denom- 
inator of the outstanding factor of the right-hand side has been 
omitted as the medium is transparent. 

Solving for R, 


A — ■ 




Ji.= ’ 


1 -~A 2 H 2 

if 


{AE v -fAHlE^ 
| AE t + i"j?irisX 




e 

m 


k*-co* 


where 
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Therefore, dropping the strokes, 


D v =icE tt -i°eAHE t , 


D. — kE, + i — eAHE, 


where 


/c = 1 + 


c 

NneA 


1 

c 2 


€ = 


NneA 


l-^-A 2 H 2 
c 1 


- = K — 1 . 


Hence, in vector notation 


D = /eE — i — eAE x H 


(47) 


Since the electric intensity is at right angles to the direction 
of propagation, V*E = 0 

and elimination of B from the field equations (1), (2), (3), (4) 


gives 


v*ve = 4d, 

a 


(48) 


the relation between D and E being specified by (47). 
Now E = E/“<**-‘>. 

d 2 E 

Therefore V* VE = — » 

dx 2 


and 


c 2 l c 2 


e , Eo + ^6^E, 


,xh} 


gia>(&c-0. 


Equating real and imaginary terms in these two equal 
expressions, 

^ = •*(5*-*^) E„ (49) 
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The second of these equations shows that if e is positive the 
plane of polarization rotates in the counter-clockwise sense when 
viewed from the source of light. Such rotation is called positive. 
Reversing the direction of the applied magnetic field reverses 
the sense of the rotation. Hence if a beam of plane polarized 
light is passed through a transparent body along the lines of 
force of an applied magnetic field, and then reflected and re- 
turned over its original path, the rotation of the plane of 
polarization is not annulled, but doubled. This magnetic rota- 
tion was discovered experimentally by Faraday in 1845. 

If the angle rotated through is denoted by a, 


da = 






. Nne* a 2 

2 jwV ; v (A 2 — a> 2 ) 2 


(51) 


where v is the index of refraction. Hence the rotation varies 
with the strength of the magnetic field and the length of the 
path. In the neighborhood of an absorption band, a> approaches 
A, and the rotation becomes very large. If the vibrating part 
of the atom is positively charged, the rotation will be positive 
when the direction of propagation is the same as that of the 
magnetic lines of force, while if the vibrating part of the atom 
is negatively charged, the rotation will be negative. Obviously, 
only the component of the field in the direction of propagation 
is effective in producing rotation. 

To find the wave slowness, eliminate E 0 from (49) and (50) 
and solve for S. Thus 

S 2 fcS 2^ 

4 tc 

or, approximately, 

s= ^(i-^£ffLy 


(62) 
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47. Metallic reflection. As a metal contains both free and 
bound electrons, the equations of the electromagnetic field take 


the form 

V-D = - G fV'Edt, 

(53) V’B = 0, 

(55) 

V X E = — i B, 

c 

(54) V x L = — (D + CE), 

(56) 

where 

D = k E, 


If 

6 = filj. 

E = E/ w(S - r -^ 


then 

V = icoS , 



d 



— = — 
dt 



and integration with respect to the time is equivalent to multi- 
plication by i/co . Therefore the field equations become 



E = 0, 



V*B = 0, 

VxL = IL+i^E. 

c\ ft)/ 


These are identical with equations (1) to (4), section 43, for 

a non-conducting medium, provided the specific inductive capacity 

is replaced by the complex quantity 

,.C 
«: + * — • 
ft) 


Therefore, remembering that the permeability is unity for 
light frequencies, the wave slowness is given by 

' + (57) 

N CD 

Hence, if S = S' 4- iS' f , 

then S ,2 -S" 2 = *cSZ, 

2S'»S" = - S 0 , 

CD 

E = E 0 e- wS// * r e lft,(s,r - t) . 


and 


(58) 
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Let the YZ plane (Fig. 22, p. 123) be the surface of separation 
between a region free from matter above this plane and a 
metallic medium below. Consider a train of plane waves inci- 
dent on the metallic surface at an angle <\. Put 

A s S 0 sin (f) t = 8 sin $ 2 , 

C ! = 3, C0S 4 * 1 ’ 

c 2 = S cos <f)g 

Then equations (40) and (45), section 45, show that 
R\\ __ j4j- 




^4 and C x being real, and C 2 complex, 
in the form 


R { 

it — np~ 

R~ pe ’ 


(59) 

This ratio may be put 

(60) 



where p and A are real. 

Suppose the incident light to be polarized in a plane making 
an angle of 45° with the plane of incidence. Such radiation 
may be considered to 
consist of two trains of 
waves, one polarized in 
the plane of incidence 
and the other at right 
angles to this plane, 
which have the same 
amplitude and are in 
phase with each other. 

After reflection these 
two trains of waves may 
have different ampli- 
tudes, and their phase 
relation may have been 

changed in such a way as to produce elliptically polarized light. 
The ratio of amplitudes after reflection is given by p in equa- 
tion (60), and the difference in phase by A. These quantities 
may be conveniently represented by means of a eTanh fFio*. 23^ 



/q 

v 

/ 

/ 

\ 

or 

\ 

or 

N 

\ 

' 7 / 

! / 

\ 



P' 

JFig. 


23 
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in which real quantities are plotted horizontally and imaginaries 
vertically. Starting from the origin O, lay off the real component 

OM=kS 2 

of S 2 y and the imaginary component 

MN = — Sg. 

a> 


From the same origin lay off the real quantities 

OQ =A\ 

QR — Cj , 

determined by the angle of incidence. Then 

QN= C?, 

and qT> = C x C„ 

is laid off by bisecting the angle RQN and making 

|3F|-Vjc?n^|. 

Connecting 0 with P and P', 

OF = A*-C l C a , 

OP =A 2 -hC\C 3 . 

Hence the ratio of amplitudes after reflection is given by 

1 0F\ 


P = 


OP I 


(61) 


and the difference in phase between the two components by 

A = ZP r OP. (62) 

48. Zeeman effect. Consider an electron which may vibrate 
under the influence of a force of restitution proportional to the 
electron’s displacement from its position of equilibrium inside 
the atom. The equation of motion of such an electron has the 


components 


d 2 x 7 
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showing that its natural vibration has a frequency defined by 


k 

~\ — 
ym 


In the presence of a magnetic field H parallel to the Z axis 
the electron under consideration is subject to an additional force 

-vxH, 

c 

so that the components of the equation of motion become 


, ell dx 
u c dt 


d 2 z 


Solving these equations, it is found that 

x = A t cos (a> t t + Sj) and y = — A ± sin (a^t + S ; ), 
or x = A 2 coa(a>/ + S 2 ) and y = A 2 sin (oajt + S a ), 

z = A s cos ( cot 4- 8), 


where 

a> 2 - — 1 // 2 = <a 2 , 


me 

or, approximately, 

+ 2 me 11 -' 

and 



me 


. e 

or 



(63) 

(64) 
( 66 ) 


( 66 ) 


(67) 


Equations (63) and (64) represent rotation in circles in the 
XY plane in the negative and positive senses respectively rela- 
tive to the Z axis. The effect of the magnetic field is merely 
to change the central force from 

kr 

to hr Refill,. 
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Consider a body which emits light in consequence of the 
vibrations of electrons which are held in the atoms by simple 
harmonic forces of the type under discussion. Suppose a mag- 
netic field to be applied in the direction of the Z axis, and let 
the source of light be viewed along the X axis. Vibrations in 
the X direction will emit no radiation in the direction from 
which the light is being observed. Vibrations in the Y direction 
will give rise to light polarized with the electric vector parallel 
to the Y axis of frequencies <o x and <» 2 , while vibrations in the 
Z direction will produce light polarized with the electric vector 
parallel to the Z axis of frequency o . ). Therefore when a source 
of light in a magnetic field is viewed in a direction at right 
angles to the lines of force, each spectral line will be resolved 
into three components. The central undisplaced component will 
be polarized with the electric vector parallel to the field, and 
the two displaced components with the electric vector at right 
angles to the field. 

If the source of light is viewed along the Z axis, no radiation 
will reach the observer due to vibrations parallel to this axis. 
Vibrations perpendicular to the Z axis will give rise to circularly 
polarized light of frequencies aq and o> 2 . Consequently, when a 
source of light in a magnetic field is viewed along the lines of 
force, each line will be resolved into two components, circularly 
polarized in opposite senses and equally displaced on either side 
of the original line. There will be no undisplaced component. 
The sense of the circular polarization of the two displaced com- 
ponents depends upon the sign of the vibrating electrons, which 
are thus shown to be negative. The ratio of charge to mass 
of the negative electron may be calculated from the displace- 
ments observed. This method was one of the earliest employed 
to obtain the numerical value of this important constant. 

While the results obtained from theory are entirely confirmed 
by experiment in many cases, a large number of lines are split 
up into more than three components by a magnetic field. It is 
believed that these are compound lines, which the optical 
apparatus employed is not powerful enough to resolve. 
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